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NOTE ON THE TRANSFERENCE OF ANGULAR MOMENTUM 
WITHIN THE GALAXY THROUGH THE AGENCY 
OF A MAGNETIC FIELD 


F. Hoyle and F. G. Ireland 
(Received 1960 March 10) 


Summary 

It is shown that, under favourable conditions, the conventional picture 
of the orientation of the magnetic field in the spiral arms of the Galaxy can 
explain the outward motion of gas at a radial speed of order 50 km sec™! 
within the central regions of the Galaxy. Such a motion cannot arise, however, 
unless at a greater distance from the galactic centre there is a corresponding 
inward motion. ‘The failure to observe such an inward motion casts doubt 
on the correctness of the Conventional picture of the magnetic field. 





1. Introduction.—_In a previous paper (1) we have considered the magnetic 
held within the galactic plane, using in our work several assumptions that have 
become conventionally accepted: 

(1) The magnetic field is orientated along the spiral arms. 

(2) The intensity in the solar neighbourhood is ~ 6 x 10~® gauss. 

(3) Initially the field had the rough form of a bar that passed through the 

nucleus. 

(4) ‘The change with time of the angular momentum of rotation in the Galaxy 

of the interstellar gas is small and can be neglected. 

Subject to these assumptions it was shown that the differential rotational 
velocity in the Galaxy winds the magnetic field into a clockspring-like structure, 
the degree of wrapping after a time interval of ~ 5 x 10° years being comparable 
with the wrapping of the spiral arms inferred by observation. 

It is proposed now to examine the fourth of the above assumptions. A 
clockspring winding of the magnetic field causes angular momentum to be 
transferred from the gas at one radial distance to the gas at another radial distance. 
Hence (4) is certainly not strictly satisfied. ‘The question is evidently whether 
or not (4) is satisfied to a sufficient approximation for the former conclusions to 
be still applicable. In the following discussion it will be shown that this 1s indeed 
the case, provided (1), (2), (3) are satisfied, except possibly at two particular 
values of the radial distance. 

2. The transference of angular momentum.—In Fig. 1 we have a spiral arm 
crossing a cylinder of radius r along the arc 4B. ‘The figure should be thought 
of as possessing unit thickness in the direction perpendicular to its own plane, 
which we take to be the plane of the Galaxy. We also take the magnetic intensity 
H to be uniform over 4B. The magnetic stress across 4B has a component 
H? cos @sin 6/47 per unit length of 4B in the direction transverse to the radius 
vector. 


19 
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It is proposed now to examine the fourth of the above assumptions. A 
clockspring winding of the magnetic field causes angular momentum to be 
transferred from the gas at one radial distance to the gas at another radial distance. 
Hence (4) is certainly not strictly satisfied. The question is evidently whether 
or not (4) is satisfied to a sufficient approximation for the former conclusions to 
be still applicable. In the following discussion it will be shown that this is indeed 
the case, provided (1), (2), (3) are satisfied, except possibly at two particular 
values of the radial distance. 

2. The transference of angular momentum.—In Fig. 1 we have a spiral arm 
crossing a cylinder of radius r along the arc AB. The figure should be thought 
of as possessing unit thickness in the direction perpendicular to its own plane, 
which we take to be the plane of the Galaxy. We also take the magnetic intensity 
H to be uniform over AB. The magnetic stress across AB has a component 


H*cos@sin 6/47 per unit length of AB in the direction transverse to the radius 
vector. 
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According to observation, the diameter d of a spiral arm varies only slowly 
withr. Inthe Galaxy d~s5ooparsecs. It will be sufficient for our purpose to take 
the arc AB to be only a moderate fraction of a quadrant. At values of r>d 
where angular momentum transfer is at all important, @ turns out to be sufficiently 
different from 7/2 for this to be the case. Inspection of photographs of extra- 
galactic systems show indeed that in many cases, even near the nucleus, AB is 
short enough to be considered approximately straight, giving AB~dsec@. The 


total force acting over AB in the transverse direction has magnitude dH? sin 6/47 
in this case. 


rs) 


‘\. SPIRAL ARM 


d \. = 


Fic. 1 


Consider the gas enclosed between two cylinders of radii r+ dr,r, the system 
again having unit thickness normal tothe plane of the Galaxy. Thenet component 
of the force acting on the gas in the direction transverse to the radius vector is 

= [(H2dsin 0), — (H?dsin8),, 4,]- (1) 
Tv 
For 6 defined as in Fig. 1, the expression (1) has positive sign when the force 
has the same sense as the rotation of the Galaxy. 

Now both d, H vary slowly withr. If 6@ also varies slowly with r, the transverse 
force (1) turns out to be small and unimportant. If @ varies rapidly with r, 
however, as it can in a certain special circumstance to be considered below, the 


force (1) can become important in promoting angular momentum transfer. In 
such a case (1) may be written to sufficient accuracy as 


-~1 wacoeo® 
—— dcosO— . dr. (2) 


The mass of gas contained between the two cylinders (again with unit extension 
normal to the galactic plane) is AB.p,.dr=dsecOp,.dr, where p, is the gas 
density. Hence the transverse acceleration is of order 


(3) 


in the sense of galactic rotation. 


Expression (3) shows the magnitude of the transverse force to be largest in 
the case =o —i.e., when the arm points in the radial direction—and when 
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|d0/dr| is as large as possible. This agrees with simple qualitative expectations. 
We consider next the special condition under which these requirements might 
arise. 

3. The rotational velocity curve.—At some initial time let P,,P, be the 
positions of two neighbouring elements of gas situated on the same magnetic 
line of force, P, at distance r from the galactic centre and P, at distance r+dr. 
If the gaseous elements move in accordance with a rotational velocity curve V(r), 
then after a time interval T we have 


ZP,00,= ©) 7, mod2z, 


r 
_ V(r+dr) 
ZP,0Q,= “your T, mod 27, 
where O is the galactic centre, and Q,, Q, are the positions after time T of the 
two gaseous elements. For sufficiently small dr, we hence have 


£0.00; = £ POP, - TS (qr (4) 


This result shows that there is no bending of the line of force at central distance 
rif d/dr (V/r) =o. This is the condition that rotation in the Galaxy shall not 
make the cos*@ factor in (3) small. We may express the condition in the form 
dV/dr = V/r, which shows that the tangent to the rotational velocity curve 
must pass through the origin. At such a point the value of @ is unchanged by 
galactic rotation. If @=o initially, in accordance with assumption (3) of the 
introductory section, then 6=o at subsequent times. (To prevent cos*@ being 
small, this assumption is not strictly necessary—all that is required is for 6 not 
to be close to 7/2 at the initial time.) 

To proceed further it is necessary to consider the likely form of the rotational 
velocity curve of the Galaxy. A plausible form is shown in Fig. 2. For central 
distance r> ~3kpc, observation indeed gives a curve closely similar to Fig. 2. 
The novelty of Fig. 2 lies in the peak at r= ~o-25kpc. Recent work by the 
Dutch radio astronomers indicates the presence of high rotary speeds near the 
galactic centre, while the work of Walker has indicated a similar phenomenon 
in M31. These high speeds presumably arise from very high central conden- 
sations of mass—a mass of ~ 2 x 10°© being required to produce the peak shown 
in Fig. 2. In this connection it is worth noting that M31 possesses an inner 
nuclear region of very high surface brightness, such as would be expected to 
correspond to a high mass concentration. A mass ~2x10°%© distributed 
uniformly throughout a sphere of radius 0-25 kpc yields a mass density comparable 
with the inner regions of globular clusters. Because of its larger dimension the 
inner galactic nucleus should therefore have a greater surface brightness than the 
globular clusters. This indeed appears to be the case for the inner nucleus of 
M3I. 

The high mass density cannot continue indefinitely with increasing r, 
however. Judging from surface brightness studies of M31, the mass density 
falls off sharply outside the inner nucleus. This implies a corresponding fall in 
the rotational velocity. Immediately outside the inner nucleus the rotational 
velocity must simply follow the Keplerian value in the gravitational field of the 
inner nucleus. That is to say, the rotational velocity falls as r-?. Sucha 
dependence is followed in Fig. 2 for 0-25 <r< ~0-75 kpc. 


19* 
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Eventually however, the field of the inner nucleus declines in importance, 
while the main field of the Galaxy increases in importance. This causes the 
rotational velocity to increase again with r. Hence it seems that the rotational 
velocity curve of the Galaxy must possess two maxima separated by a minimum. 
The depth of the minimum is uncertain, however. The value given in Fig. 2 is 
an assumption. 





i ' ' ! T T U 


ROTATIONAL VELOCITY AGAINST DISTANCE 
FROM GALACTIC CENTRE 


V (km/sec) —=— 
8 


8 








L 1 J 4 1 
4 5 6 7 8 
riiec) ——— 
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Returning now to the question of satisfying the condition d/dr (V/r) =o, 
a point of inflexion must exist on the rotational velocity curve between the 
minimum and the outer maximum. If the tangent at this point of inflexion 
intersects the abscissa at r>o, there will be two points at which d/dr (V/r)=o 
—one on either side of the inflexion. If, on the other hand, the tangent at the 
point of inflexion intersects the ordinate at V >o, there can be no such points. 
The limiting case is where the tangent at the point of inflexion passes through 
the origin. This is the case shown in Fig. 2. The condition 


s(t)- 


is then satisfied only at the inflexion. 

It is of interest to consider the values of 6 produced by the rotation curve of 
Fig. 2. This will be done for winding of magnetic lines of force over a time 
interval 7, starting from an initial state in which 8 =o everywhere along the lines 
of force. Then / P,OP,=0 in (4). The two gaseous elements remain radially 
separated by dr throughout the winding. (We consider here only the effect of 
winding, with each element maintaining its radial distance r.) Hence 


tan 0= -r7.5 (2). (s) 
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Values of 6(r) for T=5 x 10° years are given in the following table. These 
were determined using the function V(r) plotted in Fig. 2. 


r(kpc) 1°75 2 2°25-2'§ 2°75 3°0 4°0 6:0 8-0 10°0 

@ 85° 51° 0° ~— 558 81°%5 = 87°*1 = 87%4 S873 865 
Except close to the point of inflexion the field is wound in directions nearly 
transverse to the radius vector. For 2<r<3kpc we can use (3) to estimate 


the transverse acceleration, with 


rd ~ —1(kpc)- for 2<r < 2-5 kpc 


and 
dé 
a =! (kpc)! for 2-5 < r <3 kpc. 


For other values of r, @ is so close to 7/2 that our approximation used in determining 
(3)—that the spiral arm cuts the cylinder r = constant over an arc appreciably less 
than a quadrant—becomes at best only a rough approximation. But this failure 
is quite unimportant, since angular momentum transfer at such values of r turns 
out to be negligible. 

With 


d6 ~ —1 2p~ ~ —24 3 
% | 21 (kpe) » cos*@~1, p,~10-4*gmcm-, 


the transverse acceleration (3) has a magnitude of order 25H? cmsec~*, where H 
is in gauss. This may be compared with the radial acceleration, which is 


~ 4x 10!4/rcmsec™?, 


with r in cm. Putting r=10cm, i.e., r~3kpc, we see that the tangential 
acceleration is comparable with the radial acceleration if H ~ 4 x 10~° gauss. 

4. Discussion.—Subject to the assumptions (1), (2), (3), set out in the intro- 
ductory section, we can now assess the conditions required to produce the outward 
gaseous motions observed by van Woerden, Rougoor, and Oort (2). For 
2<r<3kpc it is sufficiently accurate to use the following formula for the angular 
momentum A per unit mass about the galactic centre: 


h=rV =2x 10-72, (6) 


inc.g.s. units. At r= 107 cm = 3 kpc this relation gives V ~ 2 x 10’ cmsec™, which 
agrees with the curve of Fig.2. Also we note that over the range 2<7r<3kpc we 
have V(r) approximately proportional to it, thereby yielding ho r*. 

Now for a transverse acceleration ~25H*cmsec~, acting at radial distance 
r, the change of angular momentum Ajé occurring in a time interval At sec is 
~25H*rAtcm?sec-!. From (6) we also have 


Ah~= 4x 10-7Ar, 

so that 

+25H*At~ 4x 10-Ar, (7) 
This relation requires that as r changes the gas moves always in an orbit close to 
the circular orbit corresponding to the rotational velocity curve. Since the 
transverse acceleration is applied smoothly around the orbit this should be the 
case to a tolerable approximation. The + sign in (7) must be used for 
d6/dr<o, the — sign for d0/dr>o0. To obtain |Ar/At|~s5okmsec-! we 
require H~3 x 10-5 gauss, 
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Before we go on to discuss whether or not this constitutes a plausible explana- 
tion of the results of the Dutch astronomers we note that similar effects become 
very weak indeed for r >4kpc. The differential obtained from (1) contains a 
term 


S25 Ocoee 
os d).dr (8) 


which is now considerably greater than the dé@/dr term. Using (8), the 
transverse acceleration becomes of order 


De d 
a 6 © (H2d). 
re sin 8 cos 7 (H?d) (9) 


sen Oe 
~~ > 

it can readily be estimated that (g) is numerically less than the case discussed above 
by a factor ~10~*. Hence the magnetic stress can produce radial motions only 
of ~o-5kmsec~ in the outer parts of the Galaxy. Such motions give no 
appreciable redistribu:ion of material in the time scale required to produce the 
observed degree of winding of the spiral arms of the Galaxy, viz., 5 x 10° years. 

The question now arises as to whether the form of magnetic stress considered 
above yields a satisfactory explanation of the origin of the large radial motions 
observed by the Dutch astronomers in the inner regions of the Galaxy. The 
requirement that H be of order 3 x 10~*° gauss is notasevere one. It is true that if 
the diameter d of the spiral were strictly independent of r, and if the magnetic 
flux along the arm were constant, then H at r~ 2-5 kpc would only have the same 
value as it has at the solar distance, r~8kpc, H~6x 10~* gauss. But d probably 
increases with r, implying higher values of H at smallr. Moreover, it is possible 
that magnetic lines of force escape from the spiral arms into the halo (1), and this 
would permit the flux along the arm to be greater at small r than it is at the solar 
distance. Hence it does not seem unreasonable that the magnetic intensity 
might be 5 times greater at r~2-5 kpc than it is at r= 8kpc. 

More serious is the requirement 


SIF a 
dr\rj/ ~ 


It was this condition that permitted @~o in the range 2:25<r<2-5kpc. If the 
minimum of the rotational velocity curve were appreciably shallower than it is 
in Fig. 2, the tangent at the point of inflexion would intersect the ordinate axis 
at V >o, in which case large radial motions would not arise. 

Still more serious, a large outward motion cannot arise by the present process 
without there being a corresponding large inward motion. For outward motion 
we require d6/dr <o, which occurs at r=2kpc, whereas at r= 2-75 kpc we have 
dé/dr>o. ‘Thus the motion is outwards at 2kpc and inwards at 2°75 kpc. 
This change in the sense of the mution is a consequence of the need for 
conservation of angular momentum. Subject to the assumptions (1), (2), (3) of 
the introductory section, the gaseous material at one value of r can only gain 
angular momentum at the expense of material at other values of r. This 
conservation requirement is introduced implicitly into the @-values worked out 
above, for these 6-values are only correct if sufficient gaseous material exists at 
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all values of r. The essential requirement is that to a tolerable approximation 
the gaseous material follows the rotational velocity curve. This implies that the 
gas density must be sufficiently high for the magnetic lines of force to be wound 
in accordance with the rotational velocity curve—the latter being determined 
essentially by the gravitational field of the stars. 

A radial outflow of gas in the Galaxy is observed near r=3kpc (2). Probably 
our result that the outflow should occur at smaller r, namely r~2kpc, is not an 
overriding objection to the above theory—the rotational velocity curve could 
take a form that would increase the theoretical value, whereas the observational 
value could possibly be too large. A far more serious objection is that the 
corresponding inflow of material does not seem to have been observed. It is 
true that radial velocities of recession from the Sun are observed, also of order 
50kmsec-—!, but these appear to arise from gas on the far side of the galactic nucleus 
(3) which would hence also be in outward motion. Since the observations are 
made in the 21cm line of neutral atomic hydrogen a failure to observe inward 
motion would be understandable if the inward moving material were ionized. 
There is no very obvious reason why this should be so, however. 

The failure to observe the inward motion seems indeed to be so serious as 
to make doubtful the above explanation of the outward motion. Yet within the 
terms of reference of the assumptions (1), (2), and (3) of the introductory section 
there seems to be no other way in which a rapid outward drift of gas can arise 
through the agency of the magnetic field. Hence we are led to question very 
seriously whether these assumptions are in fact correct. This question will 
form the subject of a further communication. 

Conclusions in regard to the transfer of angular momentum in the Galaxy (subject 
to (1), (2), (3) of the introduction).— The rotation curve of the Galaxy probably 
possesses at least one point of inflexion at central distance in the range 2 <r < 3 kpc. 
If the tangent at the inflexion passes close to the origin, or if it intersects the 
abscissa at r > o, there is a region between 2 and 3 kpc from the centre where a spiral 
arm remains essentially unwound. On the other hand for r<2kpc and r>3 kpc 
the winding of a spiral arm causes the angle @ made by the arm with the radius 
vector to become close to 7/2. Transfer of angular momentum is then slow 
and unimportant and assumption (4) of the introduction is satisfied. But for 
values of r in the range 2 to 3 kpc the value of | d@/dr | can be of order 1(kpc)-*. 
The transfer of angular momentum can then be very appreciable. Angular 
momentum is conserved through a compensation between inward and outward 
motion. If, however, the tangent at the point of inflexion intersects the 
ordinate axis at a value of the velocity V appreciably >o, the angle 6 is every- 
where close to 7/2, transfer of angular momentum is always small, and 
assumption (4) of the introduction is everywhere satisfied. 


St John’s College, Emmanuel College, 
Cambridge : Cambridge. 
1960 March 7. 
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NOTE ON THE BRIGHTNESS FLUCTUATION 
IN THE SOLAR GRANULATION 


3. Gaustad and M. Schwarzschild 


(Communicated by Lyman Spitzer, Jr., Director of the Princeton University Observatory) 


(Received 1960 March 21) 


Summary 


On photographs of the solar granulation obtained with a balloon-borne 
telescope, Blackwell, Dewhirst and Dollfus (1) measured a r.m.s. brightness 
fluctuation of 4:6 per cent. From this measurement they deduced a true 
r.m.s. brightness fluctuation of 18 per cent. This deduction was based on 
an accurate measurement of the contrast transmission function of their instru- 
ment, but on an approximate application of this transmission function to the 
granulation. This note aims to show that, using the same observed r.m.s. 
brightness fluctuation, the same measured instrumental transmission 
function, and the same value of 1-4 seconds of arc for the average granule 
diameter, a more detailed analysis gives a true r.m.s. brightness fluctuation 
not of 18 per cent, but of only about 7 per cent. The overestimate by 
Blackwell et al. appears to have resulted from their identifying the granule 
diameter with a wave-length rather than half a wave-length. 





Recently two rather discrepant values have been published for the true r.m.s. 
brightness fluctuation in the solar granulation, one by Blackwell et al. (1) of 18 per 
cent and the other by Schwarzschild (3) of at most 8 per cent (corresponding to 
a temperature fluctuation of 100°). The purpose of this note is to show that 
this discrepancy can be eliminated by applying a different granulation model to 
the observations by Blackwell et al. 

As a simple model we will consider the granulation as consisting of a set of 
circular disks of equal size randomly distributed over the solar surface, each 
differing from the mean brightness by a constant intensity + J). The reduction 
in the r.m.s. brightness fluctuation for this distribution will be the same as that 
for a single disk. 

The general development used here follows closely that of Uberoi (2). 


We define the two-dimensional covariance function R(é) of the brightness 
fluctuation I(r) by 


Re)=N{ [rete+ear, (1) 


where N is a normalization factor. The r.m.s. brightness fluctuation is 
proportional to the square root of R(o); 


VU) xv R(O). (2) 
The spectral density P(k) is the Fourier transform of R; 


PR)= | [ Rie * ae. (3) 
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In any telescope the image is smoothed by diffraction and other causes. 

be shown that the observed spectral density P*(k) is 
P*(k)=27S(k)P(k), 

where the sensitivity spectrum S(k) is given by 


s(k)= = [ [ ¥ Oedy, (5) 


¥O)= [ [Wow tsar (6) 


and W(r) is the instrumental profile of the telescope. (The sensitivity spectrum 
S is proportional to the square of the contrast transmission function as used 
by Blackwell et al.) The observed covariance at zero displacement is then 


R*(0)= | [ SAP CR)dk. (7) 


Finally we define a reduction factor 8 as the ratio of the observed to the true 
r.m.s. brightness fluctuation ; 


_VQP _ VRO) (8) 
V<d?> = VR) 
We now apply these general equations to our specific problem. For the 
granule model which we are considering, a disk of radius s and constant intensity 
Ij, the true spectrum P(R) is 


2 
P(k)=2nN| “Eay(hs) |, (9) 
and the true covariance at zero displacement is 
R(o) =27NI,?s". (10) 
For the instrumental profile we shall consider two cases, first the theoretical 
diffraction pattern of an ideal telescope, and second the actual measured profile 


of the instrument used by Blackwell et al. The theoretical diffraction pattern 
for an aperture of radius a at an effective wavelength A is 


W()=(ar)J,x( ) 

with 
k= —474 
ni 206, 265A 

and the sensitivity spectrum defined by equations (5) and (6) is 


k A k\?\12>)2 
S(k)= 3, =-(=)- -(:-(Z) ) | 12 
for k< ky, and S(k)=o for k>kp. 
The observations of Blackwell et al. were made with an instrument of 29cm 
diameter at an effective wavelength of 5300 A. They measured a mean granule 
diameter (2s) of 1-4 seconds of arc. By introducing these values of the para- 


meters into equations (g) and (12) we obtain the true granule spectrum and the 
ideal sensitivity spectrum. Then introducing these two functions into equation 


(sec ofarc)-', 


(tr) 
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(7) and executing the quadrature numerically gives the observed zero displace- 
ment covariance R*(o). Finally, by inserting this result as well as the value 


of R(o) from equation (10) into equation (8) we obtain for the reduction factor 
B in the ideal case 


Biaeas ial 0°76. 


In the second case the procedure is the same except that the sensitivity 
spectrum is derived not from equation (12) for the ideal diffraction pattern, but 
from the square of the contrast transmission function actually measured by 
Blackwell et al. Thus one finds for the actual reduction factor 


Bactuai = 0°63, 


which is somewhat smaller than the ideal value, but considerably larger than 
the value of 0-25 which Blackwell et al. used in correcting their measured r.m.s. 
fluctuation. 


Using this actual reduction factor, the corrected r.m.s. brightness fluctuation 
becomes 


Vv <)> =7°3 per cent, 
which corresponds to a temperature fluctuation of 
Vv (AT?) = 89°. 
This is consistent with Schwarzschild’s estimate (3) that the r.m.s. temperature 
fluctuation on the surface of the Sun probably does not exceed 100°. 

By comparing these values of the reduction factor 8 with the contrast trans- 
mission function curves given by Blackwell et al., it can be seen that the reduction 
factor is approximately equal to the contrast transmission function at a spatial 
wavelength of 2:8 seconds of arc. Thus the effective transmission function is 
approximately that corresponding to a wavelength twice the average granule 
diameter, rather than to a wavelength equal to this diameter as used in the 
approximation by Blackwell et al. 

It is clear that the result of the present note is correct only to the degree in 
which the model of randomly distributed circular disks of the size given by 
Blackwell et al. represents the actual granulation. It should be hoped that the 


present result can be superseded by a detailed analysis based on the observed 
autocorrelation function rather than on a granulation model. 


Princeton University Observatory, 
Princeton, New Jersey: 
1960 March 5. 
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Summary 


Newly determined radial velocities are given for 2 stars of spectral 
type O, 104 of type B and 14 of type A, all brighter than magnitude 7:0. 
They are part of a programme covering the region of the Scorpio—Centaurus 
moving stream, between galactic longitudes 200° and 300°. Each star has 
been classified on the revised Yerkes luminosity system. 

At least three spectrograms of each star have been measured, except in a 
few cases where previous results were merely confirmed. For 39 stars the 
statistical reliability of previously published mean velocities is increased. 
Eighteen new velocities are announced, of which nine are constant and nine 
variable. Observations have been resumed on 64 stars previously known to 
have variable velocity, including 5 double-lined binaries. 

Revised orbital elements are derived for the spectroscopic binaries 
HD 79351 and 81188, and velocity curves are presented for five separate 
cycles of the pulsating star 6 Ophiuchi. Two of these stars are members of 
the association. 





Selection of programme stars.—Blaauw (1) described in detail the motions of 
early B-type stars in the region of the southern stream, usually called the 
Scorpio-Centaurus Cluster. Earlier contributions to this study are fully reviewed 
in his monograph, and the recent literature may be traced by reference to the 
Catalogue of Star Clusters and Associations (2) and annual supplements (3). 
After establishing the stream motion, Blaauw compared the predicted and 
observed velocity components of stars with spectral types Bs or earlier. Within 
statistical limits of error, he showed that 40 stars brighter than visual magnitude 
5°5 are certain members and 50 more are probable members. For a large number 
a definite decision cannot be made until more precise values of the radial velocities 
are available. Blaauw indicated certain regions in which fainter members of 
spectral types later than Bs should be sought. 

While construction of the three-prism spectrograph for the Cassegrain focus 
of the 30-inch reflector was in progress, Dr Gollnow and Dr Buscombe planned 
the first observing programme of about 100 stars brighter than visual magnitude 6-0, 
listed by Blaauw, for which membership of the Scorpio—Centaurus stream could 
be more certainly decided if the mean radial velocity were more accurately deter- 
mined. Initial observations with the instrument during 1954 showed that 
many of these dwarf stars of spectral type B have such broad, diffuse absorption 
lines that improved velocity measures might be very difficult to obtain. In 
1955 and 1956, while Dr Gollnow’s time was largely occupied with problems of 
instrumental adjustments and photographic processing, he agreed that Buscombe 
should add some fainter stars to the programme. The wider selection included 
a number of stars of types B8 to A5 whose proper motions were found to be parallel 
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to the stream by Rasmuson (4) or Kulikowsky (5), but for which either no radial 
velocity or only a very uncertain value was available. In 1957 the entire list was 
divided so that individual observers might complete sections of the programme 
for publication as early as possible. 

Velocity system.—F ollowing Professor Struve’s personal suggestion, Gollnow 
and Buscombe selected ten certain or probable members of the cluster with reliable 
constant radial velocities to be used as internal velocity standards. These stars 
are marked with an asterisk in Table I. The mean difference in velocity in the 
sense Stromlo—General Catalogue of Stellar Radial Velocities (6) for the 81 
spectrograms is 0+0°5 (p.e.) km/s. When it became apparent that the velocity 
of + Sco (HD 149438) varies irregularly with a small amplitude, this star was 
withdrawn from further study. The internal probable error of measurement 
for a single spectrogram of a dwarf B star (usually between 3 and 5 km/s) masks 
this small variation. For giants and supergiants, which have sharper lines, the 
error is usually smaller. 

On each spectrogram of a B star between 10 and 12 lines were measured, and 
for A stars between 18 and 20. However, in the velocities of HD 105435, 120324, 
141556 and 178322 the probable errors are abnormally high, because only six 
lines could be measured on each plate. 

Ina previous communication (7) to the Society, we have outlined the technique 
of observation and measurement, our experience with the choice of absorption 
lines recommended by Petrie (8), and the general consistency of measures of 
standard-velocity stars of later spectral types. These stars are still observed and 
measured regularly, and continue to show no significant deviations from the 
adopted standard values. 

Stars with variable velocity.—For convenience in future studies of spectro- 
scopic binary or pulsating stars, the individual plate velocities for 75 stars are 
given in Table II, with references to previous work. An orbit has already been 
computed (9) for HD 170523, and observations are nearly complete for analyses 
of the velocity variations of HD 134687, 170465 and 180885. 

Table III lists new observations of the single-lined binary HD 79351 
(a Carinae), for which Curtis (10) derived a provisional orbit based on 25 spectro- 
grams obtained during 1904-07 at Santiago, Chile. The additional data make 
possible the derivation of a more accurate period, but, as Table IV shows, little 
alteration to the other orbital elements is necessary. Miss Morris has solved the 
orbit graphically by the method of Lehmann-Filhes, but the data do not warrant 
the computation of least-squares corrections. The velocity variation is shown 
graphically in Fig. 1. 

Recent observations of HD 81188 (« Velorum) indicate that the preliminary 
value of its orbital period should be increased to 116-776 days. As shown by 
Fig. 2, the fitting of both old and new measurements to the radial velocity curve 
is improved by slight changes in the longitude and epoch of periastron to 92°-60 
and JD 2416456-66 respectively. 

At the request of Professor van Hoof, Buscombe has obtained spectra of the 
B Canis Majoris-type star @ Ophiuchi (HD 157056), covering several complete 
cycles of its pulsation. The primary period is 3" 24™, with a beat period of about 
75 days (xx). Both the mean value and amplitude of the radial velocity vary from 
cycle to cycle, and a slight variation in brightness has been reported privately by 
van Hoof. References to previous spectroscopic studies of this star are listed in 
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Table II. Radial velocities from 55 spectrograms on this programme are pre- 
sented in Table V and Fig. 3. On the last night of observations a slightly distorted 
maximum of the velocity curve is seen just two cycles after the sharp one reported 
by McNamara (12) from his coudé spectrograms. It seems quite clear that 
future progress in unravelling the complex behaviour of this and similar stars will 
depend on spectra of high dispersion. We cannot report any further observation 
of Cai emission at the K line, as this region of the spectrum is not brought in 
good definition by the usual adjustment of the spectrograph. 
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Fic. 2.—The velocity curve of HD 81188. 


Further studies.—In forthcoming communications in this series, Miss Morris 
will present additional luminosity classifications for about 200 stars now being 
studied, while Dr Gollnow and Dr Przybylski will present other velocity deter- 
minations. Spectrograms are being taken of some stars fainter than 6th magnitude 
in the region of the association for spectrophotometry of absorption lines. Some 
of these objects, with proper motions and distance moduli which relate them to 
the association are being observed further for radial velocity determinations. 
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Computations are in hand for a rediscussion of the kinematics of the system as 
new, data become available. 
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Fic. 3.—Velocity curves for @ Ophiuchi on five nights in 1956-57. 


TasBLe I 


Stars with constant velocities 


Number in Henry Draper Catalogue (Harvard Annals, 91-99). 

Constellation designation. For stars lacking a Greek Bayer letter, the Gould 
number is quoted from Uranometria Argentina (Cordoba Resultados, 1). 

Visual apparent magnitude, from the Henry Draper Catalogue. 

Spectral type of the MK system, determined from Mount Stromlo spectrograms 
by A. de Vaucouleurs (V.N., 117, 449, 1957) or P. M. Morris (unpublished). 
Mean radial velocity relative to the Sun from Mount Stromlo spectrograms, to 
nearest integral kilometre per second. 

Probable error of the mean, from deviations of individual plate velocities, in 
km/sec. 

Number of plates measured. 
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Column 
8 


9 


(1) 
HD 


49336 
62747 
66591 
67536 
75149 
*77002 
*70447 
91533 
*93845 
99556 
100929 


105382 
106983 
110335 
111613 
112092 
116226 
120640 
*120709 
*120908 
*122980 


125288 
*125823 


129954 
130807 
132955 
135160 
141168 
141318 


142250 
147084 
148937 
*149438 


149711 
150898 
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TABLE I (continued) 
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Mean of individual radial velocity determinations published previously, as 
given in references. The lower case letters denote quality by internal consistency 
of the individual plate measurements, on the system used in Wilson’s Catalogue 


(6). 


The succeeding digit indicates the number of measurements included. 


The final upper case letters denote the observatories whose measurements are 
contained in the mean. 
References to publications giving previous velocity measurements. 


(2) 


Star 


30G Pup 
161G Pup 
977G Car 
79G Car 
68G Vel 
109G Car 
119G Car 
206G Car 
5?Cha 
33G Cen 
50G Cen 


92G Cen 
¢ Cru 
39G Cru 
47G Cru 
Cru 
210G Cen 
278G Cen 

3 CenA 
283G Cen 
x Cen 


v Cen 
a Cen 


22G Cir 
o Lup 
45G Lup 
40G Cir 
12G Nor 
14G Nor 


6G Sco 


o Sco 


7 Sco 


(3) 


my 


(4) 
Sp 


B3 V 
Br V 
B3 IV 
B4 Vn 
B2 Il 
B3 IV 
B3 IV 
A2I 
B3 V 
Bs IV 


(s) 
V 
km/s 
+15 
+15 
+23 

° 
+30 
+29 
+18 
+10 
+24 
+10 
+12 


+15 
+16 
+14 
—20 
+16 
a 
+10 
+ 6 


+ 7 
+9 


+10 


(6) 


p.e. 


x 


08 
2°0 
o"4 
12 
2°4 
1°4 
o's 
O'5 
2°1 
1°6 
26 


I'l 
1°2 
21 
1'°8 
2°5 
o'7 
06 
13 
°o'9 
I's 


2°0 
I°o 


28 
I'9 
18 
26 
3°0 
o'9 


o'4 
12 
1'9 
I'l 


(7) 


n 


maw AatrouwnuUw fb 


urk Pb St Qu at OOARWUMwWNHHA 


coum Nw 


(8) 


Previous 
mean vel. 


uv Qw > 


SCerTarrrenrre oreroe 


MUwWaTIAWA DLAI 


serornrororn 


nn PIAL UW SW 


one oo 


an wuw 
Fetweaete SC. see FF Peewee SFStereteee Cte re 


(9) 


References 


Lick Publ. 16 
Lick Bull. 15 
Lick Publ. 16 
Lick Publ. 16 
Lick Publ. 16 


Lick Publ. 16 
Lick Bull. 15 
Lick Bull. 15 
Mem. R.A.S. 67 


Lick Publ. 
Lick Publ. 
Lick Publ. 
Lick Publ. 
Lick Publ. 


16 
16 
16 
16 
16 


Lick Bull. 15 
Lick Publ. 16 
Lick Bull. 15 
Lick Publ. 16 
Mem. R.A.S. 68 
Lick Publ. 16 
Lick Publ. 16 
Mem. R.A.S. 68 


Lick Publ. 16 
Lick Publ. 16 


Lick Bull. 15 
Lick Bull. 15 
Mem. R.A.S. 67 


Lick Publ. 16 


Lick Publ. 16 
Ap. Ff. 64 

Ap. F. 97 & 111 
Mich. Pub. § 
Mem. R.A.S. 68 
Lick Bull. 15 
Lick Publ. 16 
Mem. R.A.S. 68 
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TAaBLeE I (continued) 


(1) (2) (3) (4) (5s) (6) (7) (8) (9) 


HD Star my Sp p.e. n Previous References 
+ 


156325 140G Sco 641 B6IV 1°3 Lick Bull. 15 
156838 46G Ara 5°88 B2vV 06 Lick Bull. 15 
Mem. R.A.S. 67 
159358 6G Ser B8 V a7 
165024 6Ara 3:90 Boss II 12 
*172910 91G Sgr 482 B3V o'8 


Lick Publ. 16 
Lick Publ. 16 
Mem. R.A.S. 68 
Lick Publ. 16 
Ap. F. 64 

Lick Publ. 16 
Mem. R.A.S. 68 
Lick. Bull. 15 


l+++ 


~ 


175191 o Sgr 214 B4sIV 28 


*175362 31GCrA- 5-41 B8&IV 1'o 


t+ | 


_ 
eopngoanages 


Oh RW COWU UW 
CWe<eweer 


Aw wnrn «WwW 


186837 73G Pav 6-42 Bs V Ir 5 


* Internal velocity standard. 
+ Kindly communicated by B. Westerlund in advance of publication. 


TaBLe II 


Stars with variable velocities 
Column 


I Henry Draper number and constellation designation. 
2 Visual magnitude and spectral type. 

3 Julian date and decimal for time of mid-exposure. 

4 Radial velocity relative to the Sun. 
5 

6 

7 


\as in Table I 


Internal probable error of deviations from the mean velocity for the plate to 
velocities for the individual lines. 

Previous mean velocity, as in Table I. 

References. 


(1) (2) (3) (4) (5) 
HD my JD V p.e. 
Star Sp ae km/s + 


57593 5°84 5086°230 +16 44 
158G CMa B3 V 6150°237 +38 Y het | 
6264-954 + 1 3°8 

6911°194 +16 2°9 


70839 6:07 6238-047 + 3 4°4 Lick Bull. 15 
87G Car B3 III 6264:991 —I1 3°4 Mem. R.A.S. 68 
6914°220 +17 3°9 
6955092 +13 1°8 


3°56 orbit: Lick Bull. 4 
see Table III 


6°06 5882°087 +9 Lick Bull. 15 
6954°141 +24 f 
6972°076 +28 
7003°976 +27 


orbit: Lick Bull. 4 
4821°979 —2I 
6236°097 +12 
6260°957 +39 
6972-095 +45 
7001°988 +38 
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TaBLe II (continued) 


(3) (4) (5) (6) (7) 
JD V p.e. Previous References 
ere km/s mo mean vel. 


4820°006 +20 4°6 +24 b 8 L_ Lick Publ. 16 
6238-128 + 4 3°2 
6265°044 +16 39 
6319°896 +12 5°3 
6321°897 +13 43 


84461 5531°065 —- 2 2°1 Lick Bull. 
154G Vel 5537°031 +13 1°4 
5826-237 +24 1°6 


84816 5527°046 5°6 Lick Bull. 
157G Vel 5919°039 2°7 
6949'210 27 
7002005 2°7 


5882°116 4°6 Lick Publ. 
5914988 4°4 

6238-146 4'I 

6263°146 

6319°908 ; 3°2 

6322°913 3°2 


88955 under observation Lick Bull. 
q Vel for orbit 


89080 4822-008 +14 ° Lick Publ. 
w Car 6322°921 ° 

6972°177 +30 

7002'036 — 2 


93237 6319°927 +49 , Lick Bull. 
24G Cha 6322°939 +42 , 
6331°921 +16 


98718 4818-094 : Lick Publ. 
7 Cen 5527°142 
6708891 


102776 4855°996 ; Lick Publ. 
j Cen 6238°198 
6263°106 
6318988 
6319°968 


5203°010 . Lick Bull. 
5527°127 
5884°201 
6238°215 
6263°118 
6319°979 
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TABLE II (continued) 
(2) (3) (4) (5) (6) (7) 


my JD V p.e. Previous References 
Sp ae km/s + mean vel. 


+23 ¢ 4L_ Lick Publ. 16 


+60 
—12 
+74 
—14 
+40 
+15 
+44 
+34 
+ 33 


° 


Porm rre 


nN oN 
ora 


AoVe . 


_ 
~ 


5938-040 


6972°204 


no > nu 
= AN CsI 0 OW 


col 


6979°154 


a 
oo 
ve 


7O002'113 


uw N 


105416 “55 6235°2 
93G Cen d ; 6263: 
6205°15 


Lick Bull. 15 


W > 


6319" 


6322: 


NNN N 


oun 


7 Ss 


105435 259° 
0 Cen 588: 


Lick Publ. 16 


+ +: 
Ow 


105937 4818: 
p Cen 6736°867 
6765°349 


Lick Publ. 16 


un Ww 


we 


107696 ° 6263°154 
24G Cru 8V 6320°004 
6323°001 

6331°946 

6736°379 


Lick Bull. 15 


“Inv NON ON 


108257 ‘980 
119G Cen 86-012 


Lick Publ. 16 
Mem. R.A.S. 67 


a 
o+ 


"955 


= 
+ 


109668 2°70 873°992 
x Mus 5265-921 
5°945 


Lick Publ. 16 


“Iw + 


111597 ‘O74 
p Cen ‘COI 
‘O14 


Lick Publ. 16 


NN NW 


112091 
p? Cru 


Lick Publ. 16 


mu 


113314 5°02 5984°980 
é! Cen Ao V 6263°192 
6265°187 
6323°038 
6356°936 
6708-988 
6730°399 


Lick Publ. 16 


NWWwN NNN 
C2O- New 


i | 
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TaB.e II (continued) 
(1) (2) (3) (4) (5) (6) (7) 
HD my JD V p.e. Previous References 
Star Sp ee km/s + mean vel. 
115823 5°53 a +30 27B + 6b 5L_ Lick Bull. 15 
202G Cen BsIII 79 -25 116A 
+Alr: : + 5 8-4 B 
Songpey +18 444A 
_ +s 53 B 
Canrene +50 r6A 
niiiini + 6 6:0 B 
33402 —16 22A 
ae ee aa 55 B 
6712002 +19 ZA 
116087 4°52 4354°083 -16 2°5 +26d 4L_ Lick Publ. 16 
208G Cen Bs V 5265°935 -— 9 6-4 
5555131 v2 3°9 
5995°947 — 5°° 
118978 5°53 5893°216 +21 2°3 -—30c 2L_ Lick Bull. 15 
251G Cen Bg LV 5994°949 —20 2°9 
6319°069 - 4 2° 
6323°051 —- 2 2°4 
6718°005 + 8 1°4 
6736°951 —14 3°3 
120324 3Z°2V 5241°004 — 5 4°6 +13 b 5 L_ Lick Publ. 16 
pe Cen B3 Ve 5268-921 — 4 Cig. 
5588'035 +II 8 13°2 
5995°963 —14 5°0 
121190 5°84 6263:230 +31 52 + 8b 4L_ Lick Bull. 15 
288G Cen B38 V 6265:201 +38 58 
6320°032 +42 50 
6323°071 +37 3°9 
121743 3°86 4827°165 +11 21 +7b 7L_ Lick Publ. 16 
¢ Cen B2 V 6004°940 — 16 3°3 
6712°037 +19 5‘! 
6766-901 — II 7°3 
121790 3°92 5266-928 +14 370 + 7c 7L_ Lick Publ. 16 
v' Cen B3 IV 5268-927 — 3 3°0 
5586-088 + § 4'0 
5650°938 + £ 5‘0 
6649°209 +48 7°6 
6712°053 — #3 5'1 
6766-864 + 9 5°3 
6766-905 — 44 6-0 
123335 6°43 5265°956 —14 40 + 3c 3L_ Lick Bull. 15 
318G Cen Bs IV 5585°084 +39 5°4 -18 c 3 R Mem. R.A.S. 67 
124367 5°20 5203°124 - 1 64 +19 d 6L_ Lick Publ. 16 
328G Cen B3 Ve 5290°94I + 2 6°1 +11 d 4R Mem. R.A.S. 67 
5588-071 +21 6°5 
5998-968 88 4°7 


20* 
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TaBLe II (continued) 


(3) (4) (5) (6) (7) 
JD V p.e. Previous References 
er km/s + mean vel. 


§265°982 — 2 5d 5 Lick Publ. 16 
5591°128 +28 . 4¢ @ Mem. R.A.S. 68 
5995°990 +18 : 


4826-201 ‘ Lick Bull. § 
5§260°112 ‘ Lick Publ. 16 
7002°217 . M.N. 116 


4858-124 . Lick Publ. 16 
6355°995 
6709°059 
6717°994 
6765:°868 


131492 5267-001 ° d Lick Publ. 16 
@ Cir 5290°954 , Mem. R.A.S. 67 
5586-129 
5998-981 


132058 . 4822°227 ° Lick Publ. 
6 Lup 5268-965 
5585°139 
6711°023 


133242-3 5267-017 ° Lick Publ. 16 


a Lup 5268-999 
5591°140 + 1 
5996°033 ° 


133955 5267-028 +17 . Lick Publ. 
A Lup 5269-003 + 8 , M.N. 119 
5591°167 +? 
5998-992 oy 


134687 under observation : Lick Publ. 
66G Lup for orbit M.N. 116 


135737 6356-019 +21 
7G TrA 6736-984 —13 
6738-967 —12 
6798-876 —10 


136504 se -— 9 ‘8 E Lick Publ. 
« Lup 5297°964 428 


‘Qa —19 
6793°876 ae 
6798-856 iy 
137432 ° 5269°014 +18 , Lick Bull. 
100G Lup 5651-008 8 


3$999°014 
6322°137 


138690 : 4971°879 . Lick Publ. 
y Lup 5996-058 M.N. 119 
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TaBLe II (continued) 


(1) (3) (4) (5) (6) (7) 
HD JD V p.e. Previous References 
Star 5 ee km/s + mean vel. 


b Lick Publ. 16 
c Pub. D.A.O. § 
a Ap. f. 64 

b Ap. F. 42 


138764 521r°m51 + 4 5 
90G Lib 5999°030 —12 3°7 8 
: — 64 3.2 2 
6720°020 { +13 od 4 
6736006 —- 5 4°0 
748-960 ° 3°5 
139160 6:03 6356°997 +26 92 
94G Lib Bg IV 6392°947 +28 8-5 
6718°043 + 4 4°7 
6765°894 —12 6°3 


139365 3°63 4829°215 + 5 77 L_ Lick Publ. 16 
7 Lib B4 V 5586-160 +44 5°2 ; Ap. F. 64 


141556 4°11 4826-251 + 8 2°6 Lick Publ. 16 
Lu Ao III , 

— _Ty _-5555°246 mid 

2°2 

5613°086 “ 

5897°285 sa 


141637 4°65 5211°161 52 Lick Publ. 16 
1 Sco B3 V 5612-092 2°9 
5996-087 7:0 
6263°272 8-1 
142139 5°76 6354069 2°6 
19G TrA A1V 6720°046 19 
6737°023 19 
739°004 2°0 
6765°923 2°3 

142165 5211°177 4°0 Lick Publ. 16 
4G Sco 6005°999 5°4 
JII‘III 4°6 
6736-018 3°4 
748-972 4°4 


5586-171 58 Lick Publ. 16 
5612°139 30 
6006°026 6°9 
6736°030 71 
6766-966 5°6 


142990 6319°174 4°9 Lick Publ. 16 
11G Sco 6323°152 37 
6357°049 48 
6768-929 72 


145482 . 5635°048 50 Lick Publ. 16 
13 Sco 5646°047 3°9 Ap. Ff. 115 
5996°099 4°2 
147971 . 4971911 3°3 Lick Bull. 9 
e! Nor 5291°032 4°2 
5313°938 5°3 
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TaBLe II (continued) 
(3) (4 (5) (6) (7) 


JD V p.e. Previous References 
EN km/s + mean vel. 


6037°954 - -19b 4L_ Lick Publ. 16 
6320°131 
6322°189 
6326°154 


Lick Publ. 16 
Ap. F. 7o 


148605 5211°201 
i Sco 5555°290 
6003°073 

6035°964 

6319°188 

6323°173 


i ee eae eG 
mw ons ad 


148703 "25 §291°054 Lick Publ. 16 
72G Sco F 5292'996 
6709°124 


6737°044 


Np Ww 
Owe NN 


149757 2°56 4943°967 
COph Oo9.5 Vnk = 5296-025 
5647°038 


Ap. F. 64 
Lick Publ. 16 
Pub. D.A.O. § 


oO AN 
ws 


153613 5°06 5623°109 Lick Publ. 16 
117G Sco B8 V 5914°280 
5995°091 
6379°042 


wun 


wm w 
NW O 


o 2) 


154481 6°20 5350°909 
80G Oph AolIII-IV 6037-972 
6720°089 
6737°066 


~ NH & N 
a3 = tn 


155450 6:00 g ‘219 3 Lick Bull. 15 
130G Sco Br III 5°01 + 37 2 Mem. R.A.S. 67 


5°308 


u 
° 


fe) 


127 


* 
fo) 


157056 x see Table V s Pub. D.O. 8 
@ Oph j Ap. 7. 64 
} Lick Publ. i6 
Ap. F. 124 & 128 
P.A.S.P. 69 


157243 Lg 3°13 : F@e 4 Lick Publ. 16 
146G Sco 


6711°17 


157246 . 4848-240 ; : Lick Publ. 16 
y AraA 4943°977 
5266-062 
5291-010 
5647°049 
5999°075 
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TABLE II (continued) 


(3) (4) (5) (6) (7) 
JD V p.e. Previous References 
ee < * mean vel, 


5914°309 39 +18 c 8L_ Lick Publ. 16 
5919°308 30 
6709°171 Z°1 
6737°084 2°3 


4965°935 72 Lowell Bull. 1 
5293°052 6-2 Lick Publ. 16 


4944'001 3°4 Lick Publ. 16 
4970°913 2°9 
4973912 21 
5269089 43 


6036-021 373 
6037°993 6-2 
6320°183 4°2 
6322°224 
6720°125 


5586°257 Lick Bull. 15 
5995°144 + 34 ' 

67138:°160 +39 

6737°100 a 

6768-051 — 9 


5627°160 +29 
6379°090 —10 
6720°154 + 6 
6804°940 —- 7 


under observation Lick Publ. 16 
for orbit 


Lick Publ. 16 
Mem. R.A.S. 67 
orbit: M.N. 116 


5290°078 +10 7 Lick Publ. 16 
77G Sgr 6323°208 +29 

6353°160 +14 

6356°180 +38 

6393°077 + 9 


173117 . 5612°212 
94G Ser 5985°171 
5987°173 

§997°167 


173948 , 4856-286 , +20 b 4L_ Lick Publ. 16 
A Pav 4966°057 . +14c¢ 4R Mem. R.A.S. 68 
5293°108 , 
5314°006 
5999°117 
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TaBLe II (continued) 


(3) (4) (5) (6) (7) 
JD V p.e. Previous 
SOR ceux km/s + mean vel. 


4906-026 — 3 , +6d 3L_ Lick Bull. 15 
5260°172 — 38 ; 

5314°026 +30 

§715°947 = 

5996°176 —I5§ 


References 


180885 ° under observation Lick Bull. 15 
162G Sgr for orbit Mem. R.A.S. 67 


182180 5999°153 +14 
179G Sgr 6036°071 — 3 
6038-046 +30 

6066-988 —19 


91°133 +s ‘ Lick Publ. 16 
56-948 “el 
37°181 +15 
68-087 —10 


2 
3 
7 
7 


5 
5 
6 
6 


Tase III 
Radial velocities of HD 79351 


JD V 
ee km/s 


4793°058 + 6 
5585°890 +45 
5826-216 +20 
5882°052 +30 
5894°030 +18 


6238-081 +20 
6263°040 + 9 
6318-938 +26 
6331°903 +20 
6647°969 +3 


TaBLe IV 
Orbital elements of HD 79351 


Provisional Revised 
69-744 64-74469 
+233 +23°3 km/sec 
21°5 21°5 km/sec 
JD 2416533°81 JD 2416534°215 
o'18 o'18 
115°°84 112°:97 
a sini 1°96 x 10° 1°96 x 10° km 
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TABLE V 
Radial velocities of 6 Ophiuchi 


p.e. JD 
+ i 
2°6 5983°085 
"099 
113 
5613140 ‘ "126 
"142 
5625°071 : ‘155 
085 } ' *168 
"100 ‘ 182, 
113 
"126 
5987-099 
"140 ‘141 
“154 
"168 
“182 5994-981 
*196 ° 5995°006 
‘210 . 034 
054 
076 
5645°933 
‘960 . 098 
‘982 ' 133 
5646-011 ° "154 
‘O31 


"055 : 6034'913 

085 ; "931 

"102 054 
128 : ‘973 
‘989 

5983°016 . 6035°015 
"030 : 032 
043 . 056 

058 - -067 

‘072 . 083 


+++ it ++ti 
wuwnne NM HON HN 


* Under-exposed 


TaB.Le VI 


Observers and measurers 


Name Observed Measured 
W. Buscombe 310 411 
H. Gollnow 88 ° 
G. Hagemann 2 42 
W. Heintz* 16 ° 
P. Morris 68 
A. Przybylski 105 


Totals 


* Munich Observatory: guest observer, 1954-55. 
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THE EMISSION SPECTRUM OF THE NIGHT SIDE OF VENUS 
B. Warner 


(Communicated by the Director of the University of London Observatory) 
(Received 1960 April 5) 


Summary 
The measurements given by Kozyrev of features in the spectrum of the 
night side of Venus are examined. The existence of nitrogen bands is 
confirmed and evidence is presented identifying emission lines of neutral 
and ionized oxygen. 





In a recent paper (1) concerned with the night side of Venus, N. A. Kozyrev 
lists 50 features observed by him in the emission spectrum of the Ashen Light. 
Of these, 22 features are tentatively identified as being produced by the N, and N,*+ 
molecules, the others being left unidentified. A further study of the Kozyrev 
spectrum has shown that many of the latter may be due to atomic oxygen emission. 

1. Kozyrev’s spectrum.—The spectrum obtained by Kozyrev contains features 
that are at the limit of detection of his instruments; consequently the reality of 
many of these may be questioned until further evidence is produced. An attempt 
to confirm the Kozyrev spectrum has been made by Newkirk (4) and is briefly 
described in Section 4. In the present paper it is necessary to accept the reality 
of the whole spectrum as Kozyrev gives no indications of the objectivity of 
individual features. 

In Kozyrev’s paper, a list of wavelengths of spectral features is given which 
distinguishes between ‘lines’ and ‘bands’. In fact, it is evident from the 
reproduction of the spectrogram that the ‘lines’ have a width of several angstroms 
and there is a gradual transition in width from ‘lines’ to ‘bands’. The latter are 
better termed emission ‘strips’ since the features identified by him (N, and 
N,* band heads) coincide with ‘lines’ in his spectrum. 

From the data given in (1) the line-width due to the slit is found to be 1-25 A, 
while the photographic resolution could not be less than about 2-5 A. From the 
reproduction and data it would appear that the actual resolution is 4-5 A, thus 
suggesting that the emission strips may be due to unresolved lines and hence can 
be caused by atomic transitions. 

Accordingly, a search has been made for atoms capable of producing the 
observed features. To do this, Kozyrev’s measurements are compared with 
wavelengths taken from Moore’s tables (2). ‘To test whether the number of 
coincidences between these two tabulations is significant the number expected on 
chance is computed as shown in Section 2. For completeness the coincidences 
with bands in the N, and N,* spectra are also treated quantitatively, taking data 
from Pearse and Gaydon (3). ‘Throughout this investigation only strong lines 
and bands in the listed spectra are considered. 

2. Statistical considerations.—We shall count a coincidence whenever a labora- 
tory line falls within x units of an observed line. ‘Then, if there are M observed 
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lines within a spectral region of X units, the chance of coincidences for one laboratory 
line from the same region is 
2Mx 
X 
If there are N lines in the laboratory spectrum then the total number of 
expected coincidences C is given by 
C=Np 
2MNx 
: a 
For emission strips the corresponding formula is 
2M'N’ (x+ 3w) 
X’ 
where w is the average width of the strips (10 A) and x is measured from the edge of 
a strip (A). The average interval between lines (X//) is 44 A and the average 
interval between strips (X’/M’) is 61 A. 

In Tables I to V an asterisk (*) is placed whenever the number of observed 
coincidences (O, O’) is significantly greater than C, C’ ata 5 percentlevel. Fora 
real correlation, O, O’ should remain nearly constant as x increases. From a 
study of the results given in Tables I to V this is found to be the case for O in Tables 
I and II and O’ in Tables IV and V. 

3- Comparison of numbers of observed and calculated coincidences 

(i) Nitrogen.—The three important systems of bands in the N, and N,t 
spectra are the Main System (MS), the Second Positive Group (2PG), and the 
Vegard-Kaplan System (VK). These are treated separately and compared first 
with lines (O and C) and then with emission strips (O’ and C’) in the Ashen Light 
emission spectrum. ‘Tables I to III give the results. 


C= 





TABLE I TaBLe II 
N, (2PG) N,* (MS) 
ty oo € 
o'6 * 06 
13 - a°3 
I‘9 1°9 


2°5 2°5 


~ 
~ 


wnNnNN O 
ww NNN O 
us oO O 


TasLe III 
N, (VK) 
“a c 
- 18 
"4 2'2 
o'8 25 
1°2 29 
:. 26 3°2 
It is clear that the coincidence of N, (2PG) and N,+ (MS) bands with observed 
lines is significantly above chance and O is approximately constant as expected. 
There is no evidence for the emission of the N, (VK) system. Both of these 


results confirm Kozyrev’s results. The N, and N,* bands do not seem to coincide 
with the strips. 
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No other molecules have been found that give significant results, so the search 
has been extended to atomic lines with the justification of Section 1. 

The only atom found that gives any significant coincidence is oxygen. 

(ii) Oxygen.—The 4 O1 lines of intensities 8 or greater and the 31 O11 lines of 
intensities 7 or greater (within the required spectral region) have been compared 
with Kozyrev’s list. 'The numbers of coincidences are given in Tables IV and V. 


TaBLe IV TABLE V 

O!1 Ou 

’ Cc oO’ 
= 12* 
14 13* 
2°8 3° 
4°2 13* 
5°6 13 
70 14 10:2 
8-4 14 I1°2 


Table V is more extensive than the previous ones as the evidence for the pre- 
sence of oxygen rests mainly on it. 

It appears from these two tables that there is some evidence for the presence of 
O1 and strong evidence for O11, the latter contributing mainly to the emission 
strips. The O’ for O11 are significantly greater than the C’ for small x and O’ is 
virtually constant even to x=6A. 


TasBLe VI 


Identification of features in the emission spectrum of Venus 


Venus Identification 
Lines and 
strips Probable Possible Absent 
3528 3532°6(MS)([1,0] 
3536 3536°7(2PG)[8] 
3548-9(MS)[3, 0] 
3563°9(MS)[4, 0] 
3577 
3583 
3616-3630 
3646 
3655 
3665 


3576°9(2PG)[10] 
3582'1(MS)[4, 3] 


3641°7(2PG)[3] 


N- = =» NN 


3671°9(2PG)[6] 
3686 
3700 
3710-3716 
3724-3732 
3746-3756 
3781 
3805 
3825 
3836 
3860 


3712°8 O 11[7], 3710°5(2PG)[8]} 
3727°3 O1(8] 
3749°5 O19], 3755°4(2PG)[10] 


3804°9(2PG)[10] 
3823°5 O1fro] 
3835°4(MS)[1, 1] 


Ns & & NNN HN 


3857°9(2PG)[s5] 


3888 3884-3(MS)([3, 1] 


3882:2 Onf[7] 


3894-6(2PG)[7] 
3900 
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TaBLe VI—continued 
Venus Identification 
Lines and 
strips 
3914 3914°4(MS)[6, 0], 3912°0 O1l[r1o] 
3946 3 394773 O1[10], 3943°0(2PG)[3] 


Probable Possible Absent 


3954°4 O11[7] 


3958 39547 O1[io] 
3973°3 Ouro] 
4000 3998°4(2PG)[9] 

4024-4035 

4052-4056 4059°4(2PG)[5] 


_ 72:2 O1i[8] 


075°9 O1i[10] 


4094°8(2PG)[4] 
4105'0 O11[7] 
41192 On[8] 


‘8(2PG)[5] 4140(MS)[2, 0 
4153°3 O1i{7] 
4166°8(MS)[3, o] 
4184-4192 85°5 O18], 4189°8 O1i[10]} 
4199°1(MS)[4, 2] 
4200°5(2PG)[6] 
4200-4219 
4236°5(MS)[7, 5] 
4250 3. 4253°7 O1i[S] 
4264 
4275 
4314-4330 
4340-4350 
4363-4370 
4390-4395 
4410-4420 
4430-4440 
4450-4462 
4496-4504 


4269°7(2PG)[5] 


Ww & 


78:1(MS)[8, 6] 
431771 O11[8], 4319°6 O11[5] 
4345°6 O1[7], 4349°4 O1[5] 
4366-9 O11[7], 4368-3 O1[10] 
4396°0 O1[7] 
4415°0 O1r{10], 4417-0 Orn[S8] 4416°7(2PG)[3] 


we 


NR N NN Nw WN WN 


4515°9(MS)[6, o] 
4554°1(MS)[4, 0] 
4591°0 O11fg] 
4596-2 O11[$] 
4599°7(MS)[6, 2] 
4012-4616 

4640 ‘8 Onf[g] 

4652 51°8(MS)[4, 13], 464971 O1t[10] 


4670 4676-2 O1[8] 


4661°6 O1[g] 


4699°2 O11[7] 

4705°4 O11[8] 

4710-4720 2 4709°'2(MS)[4, 2 
4760-4780 2 

4856 2 

Having established, on statistical grounds, the probable presence of certain 

nitrogen and oxygen emission features it is possible to construct a list of identifica- 

tions—Table VI. In this table the first column is a reproduction of Kozyrev’s 


published wavelength measurements, which have a probable error of 2-3 A. 
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Single numbers denote ‘lines’, and pairs of numbers, ‘strips’. The second 
column gives his estimates of the intensities of the features. Column three is a 
list of probable identifications. The square brackets signify the intensities—for 
the M.S. the first number gives the intensity in a discharge through helium con- 
taining a trace of nitrogen, the second the intensity through nitrogen at low pressure 
(3). The oxygen and nitrogen intensities are not on the same scale. The next 
column gives possible identifications and the final column absent features. 

The reason why oxygen and not nitrogen causes apparent ‘bands’ in the 
emission spectrum is due to the resolving power of the optical system used and the 
multiplet structure of the oxygen spectrum. It can be seen from Table VI that 
all but one of the pairs of features that cause emission strips have separations less 
than 4°5 A, i.e. less than the probable resolving power of Kozyrev’s instrument. 

4. Conclusion.—By statistical means we have found evidence for the pre- 
sence of Ns, N,* and O11 in the atmosphere of Venus, and shown the presence of 
O1tobe probable. The absent strong lines (see last column of Table V1) are a 
cause for concern. However it should be remembered that the detection of all 
the features in the Ashen Light spectrum has been marginal and some detection 
failures must be expected. Moreover, differences of relative intensities in the 
laboratory and the Ashen Light might result in the non-appearance of some lines 
whose laboratory intensities are thought to be above threshold. In the Ou 
spectrum all but three of the lower excitation lines (low E.P.<24eV) and only 
four of the higher excitation lines have been detected. This suggests that the 
laboratory excitation temperature was higher than that of the Ashen Light. If the 
lower excitation O11 lines only had been used in the statistics the calculated 
significance of O11 presence would have been much greater. 

The attempt to confirm Kozyrev’s work by Newkirk (4) provides no additional 
evidence in favour of oxygen, the one extra feature discovered by him still being 
unidentifiable. Newkirk provided a criterion for deciding the reality of the 
emission features in his spectra, but this throws no light on the fainter features in 
Kozyrev’s spectrum as Newkirk’s apparatus had a lower dispersion. Two of the 
brighter bands (due to N,*) are confirmed. 

One might conclude that the reality of the emission features in Kozyrev’s 
spectrum is much strengthened by their independent identification with oxygen 
lines. 


The writer’s thanks are due to Mr T. Kiang for advice relating to the statistical 
problems concerned in this paper. 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7: 


1960 April 4. 
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THE DETERMINATION OF THE INCIDENT FLUX 
OF RADIO-METEORS 


T. R. Kaiser 


(Received 1960 April 7) 


Summary 


The observed rate of radio-meteors depends on a number of factors 
which include: (i) the incident flux and distribution in magnitude of the 
meteors, (ii) the radiant co-ordinates, (iii) the variation in ionization along a 
meteor trail, (iv) the nature of the radio reflection process and (v) the 
parameters of the radio-echo equipment. ‘The present work represents a con- 
siderable simplification over previous attempts (5, 6) to relate the radio echo 
rate to the actual flux of meteoroids. Simple formulae are derived relating 
these two quantities which, when used together with suitable graphical 
representation of the directional properties of the aerial system and of the 
reflection geometry, enable the incident flux of shower meteors to be deduced 
from the observed rate. The fact that the variation with time of the shower 
rate obtained with a fixed aerial beam may be predicted is relevant in con- 
nection with a method (7, 8) for shower radiant determination. Limitations 
due to the approximations in the theory, and to the simplified ablation theory 
on which it is based, are discussed. 





1. Introduction.—The study of radio reflections from meteor trails has 
yielded a considerable amount of information concerning the properties and 
distribution of interplanetary matter and the physical processes accompanying 
the entry into the atmosphere of meteoric particles (meteoroids). A virtue of 
the technique is its inherent objectivity; the mechanism and geometry of 
reflection are fairly well understood and the essential equipment parameters 
(transmitted power, receiver sensitivity, aerial gain, etc.) are amenable to precise 
determination. It may therefore seem at first sight surprising that relatively 
little effort has been devoted to the problem of deducing the absolute incident 
flux of meteoroids from the observed echo rate. The following are the major 
difficulties which arise: (1) The reflection process is selective; echoes are only 
observed from trails which produce sufficient ionization in the vicinity of the 
specular reflection point (where a vector from the observer makes a normal 
intersection with the trail). (2) The sensitivity varies with direction in the 
aerial beam and with range. (3) In the case of sporadic meteors the relation 
between incident flux and echo rate depends on the orbital distribution. In 
view of the additional complication resulting from (3) the present discussion 
is limited to the case of shower meteors. 

In former theoretical work (1,2) concerned with the distribution of 
reflection points it was shown that they should be confined to a relatively narrow 
height interval about a mean height near 100 km and in spite of certain simplifying 
assumptions (particularly as concerns the meteor ablation process) fair agreement 
with observation has been obtained (3, 4). A previous extension of the theory 
to the present problem (5, 6) suffered through involving unduly complicated 
integrals which only in limiting cases (e.g. for a narrow, gaussian aerial beam) 
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led to manageable formulae. In the present work simple formulae are derived 
which combined with a suitable graphical procedure enable the magnitude 
distribution and incident flux of shower meteors to be derived from the observed 
echo rates. The theory enables the variation in echo rate due to the diurnal 
motion of the meteor radiant to be predicted and hence it has some importance 
in connection with the method for shower radiant determination developed. by 
Clegg (7), particularly as modified by Keay (8) for high echo rates. 

2. The geometry of meteor reflections—On account of the specular condition 
all reflection points lie in the ‘echo plane’ ABCD (Fig. 1), which is normal to 
the radiant direction and passes through the observing station, O. The elevation 
of ABCD above the horizontal at O is thus equal to the radiant zenith distance, 
Xo» measured at O. The reflection points are further limited to a band 
corresponding to a relatively narrow height range about a mean height 4 which 
defines the locus SS! in the figure (we will call SS! the ‘echo line’). 


Zenith 


lant 











D 
Fic. t.—The echo plane ABCD which is normal to the radiant direction and contains the echo 


line SS'. 


The minimum detectable line density, «,, in a meteor trail which intersects 
the echo plane at P is related to the system parameters by (5, 9, 10)*. 


4272 3 12 2 
tp = (Far) (= ) for a»<2-4x 101? cm-! (decay type echoes) (1) 


PG e 


and 
= 4 3p 
tp = (ae) (=) for a>24x10"%cm-! (persistent echoes). (2) 
PyGa e 

Here, R=OP, P,=minimun, detectable echo power, P,=transmitted power, 
G=V\ G,G, where G, and G, -re the transmitting and receiving aerial gains 
respectively, A= wave-length, m, e, and c are respectively the mass and charge 
of the electron and the velocity « f light (in c.g.s. gaussian units). 

From (1) and (2) we may deduce the variation in sensitivity for meteor 
detection, i.e. of ap, along the echo line SS'. However, since the echo plane 
moves through the aerial beam due to the daily motion of the radiant, it is more 
convenient to use the representation of Clegg (7) in which the contours of 
constant «, are plotted on a sphere (the ‘echo surface’) at height 4 above the 
Earth’s surface. ‘This is illustrated in Fig. 2 for one of the rhombic aerials at 


* The echo enhancement due to plasma resonance (g) with transverse polarization has been 
ignored in(1). ‘The numerical factor in (2) has been modified in accordance with a correction to the 
original theory due to Manning (11). 


2! 
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Sheffield with h=100km. The projection is such that the radial coordinate 
is proportional to the distance R, measured in the echo surface from its inter- 
section with the observer’s zenith; the normalized sensitivity contours are 
represented by lines of constant p=apo/%p, where apo is the minimum value 
of a» and depends on the equipment parameters as well as the aerial gain*. 
The error in scaling distances in the curved echo surface directly from the 
projection of Fig. 2 is negligible. 


(ee, 


500 —e Rgkm 1000 
200 yO 400 500 600 700 400 900 00 
—> Rkm 
Fic. 2.—Normalized meteor sensitivity contours for one of the Sheffield rhombic aerials, pro- 
jected onto the echo surface at 100 km height. Rg=radial distance measured in the echo surface 


from its intersection with the observer's zenith, R=slant range, ¢,=elevation, p=azimuth, 
Ga = azimuth of beam axis (azimuths measured E of N). 
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Fic. 3.—Echo lines on the 100 km echo surface (similar projection to Fig. 2). 


The echo lines may also be represented on the echo surface as shown in 
Fig. 3. Given the radiant zenith distance (xy)) and azimuth (gr), Fig 3 can be 
superimposed on Fig. 2, and oriented so as to dispose correctly the radiant 
azimuth relative to that of the beam maximum, thus giving the intersection of the 
echo line with the sensitivity contours. In this way the variation of a, along an 


* Provided apy < 10"? it is sufficient only to employ (1) in evaluating p, since in this case the con- 
tribution of persistent echoes to the total rate will be small. 
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echo line is obtained. The Geminid echo lines at various times relative to radiant 
transit are shown in Fig. 4 for an observer at Sheffield (latitude 53° 23')*. 

If the echo rate is sufficiently high, considerable simplification may be achieved 
by considering only echoes within a restricted range interval centred on some 
range R (e.g. the range at which p=1). In this case a set of curves of the type 
shown in Fig. 5 will be useful. These are loci of constant azimuth and zenith 





Fic. 4.—Geminid echo lines for the latitude of Sheffield (53°23'N) in a similar projection to 
Figs. 2 and 3. The figures on the curve are (T, xo) where T hr. is the radiant hour angle and xo 
degrees is the zemith distance. 


90° 























Fic. 5.—T=hour angle at which the echo line for a radiant at declination 8 intersects azimuth 
WE of N) at slant range R=450 km ($,=10°'9°, latitude 53°23'N). The broken curves are loci 


of constant zenith distance, xo. 


distance of a point on the echo line at the fixed range R, plotted in radiant 
declination—hour angle coordinates. Given a set of such curves, together with 
the sensitivity contours of Fig. 2, the sensitivity factor p and zenith distance x, 


* Formulae relating xo and ¢y to the celestial coordinates of the radiant and for calculating the 
echo lines in the projection of Fig. 3 are given in the Appendix. 


21* 
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can be determined as a function of hour angle for any specified radiant 
declination. In particular the time at which the echo line intersects the point 
of maximum sensitivity (on the beam axis) is easily ascertained. A method for 
computing these curves will be found in the Appendix. 

3. The representation of the incident flux of shower meteors.—An approximate 
relation for the electron line density «cm~-! produced by a spherical 
meteoroid travelling through a rarefied isothermal atmosphere is (12): 


“= am(p/Pm)(1 —}p/Pm)? (3) 
where the maximum line density «m occurs at pressure Pm; %m is proportional 
to the initial mass of the meteoroid and to cosy. If the pressures p, pm occur 
at heights h, hm, then 

h—hm= Hn (p/pm) (4) 
where H is the exponential scale height. (3) and (4) define the line density 
as a function of height (Fig. 6). Weiss (13) has obtained a more exact equation 
and has considered the effect of a linear height gradient of H. However in neither 
case is the form of the electron density profile altered to an extent which would 
be serious for the present purposes. Indeed major departures from the simple 
ablation theory are more likely to come from processes such as fragmentation (14). 


- 
| 


Bite. ee 


i i 4 
<a n-H . 
H 
Fic. 6.—Meteor ionization profile. 





In view of the dependence of am on x it is convenient to introduce the 
quantities 
az = &m (cos y)~? (5) 
which we will call the maximum ‘“‘zenithal’’ line density (i.e. the line density 
which the meteoroid would produce if incident vertically) and the absolute 
radio magnitude M which is defined by Kaiser (10) as 
M = — 2:5 logio%z +35. (6) 
Decay type echoes occur for M> +5 
The incident flux of meteoroids (per unit time across unit area normal to 
the radiant direction) can thus be represented by: 
(az) da, = flux producing zenithal line densities between «, and a, + daz, 
©(«,) = flux producing zenithal line densities greater than xz, 
A(M) dM = flux producing meteors with magnitude between M and M+ dM, 
B(M) = flux producing meteors brighter than magnitude M. 





No. 3, 1960 The incident flux of radio-meteors 289 


In many cases we can approximate to the distribution function ® by the power 
law 


(az) =caz-8 (7) 


where c and s are constants. Such a law appears to hold for sporadic meteors 
over a wide range of magnitudes with s close to 2:0 (2,5). In a study of the 
major showers (6) s was found to have some dependence on «, and to be 
generally less than 2, at least for the sub-visual radio meteors. If (7) is valid 
we obtain: 

c 


O(%2)= (8) 


and 
B(M)«A(M)oca™ (9) 


where s=1+2-5loga. For s=2 we get a=2°51 which corresponds to a constant 
total mass of meteoric material per unit magnitude interval. 

4. The echo rate per unit length of echo line.—With the help of the geometrical 
considerations in Section 2 we can now relate the observed shower echo rate to 
the actual incident flux of meteoroids. 

Consider those meteoroids incident upon a strip of the echo plane at P (Fig. 1) 
which is normal to the echo line and is of width dl, and let the height and 
atmospheric pressure at burn-out be hy and py respectively. Equation (3) may 
thus be re-written: 


a/am = raat —u) (10) 


where u = p/py=exp (—x/H),x=h—hy. The number of meteors passing through 
the interval of the strip between heights x and x + dx above /, and with maximum 
line densities in excess of «m is simply 


an =0( am \ oa (11) 


cos x/ sinx- 


The radio-echoes from these will be detected provided the line density at the 
echo plane is in excess of ap, i.e. provided* 


t -1 -2 
> ome a 


The total number of echoes obtained from meteors intersecting the strip is 
therefore 
1 
N,dl= 2 (*9¢,.)%. (13) 


sinx Jo u 


Where «,,= —?- u-1(1—u)-*; it is the maximum zenithal line density 
27 COs x 
of the faintest detectable meteor at height hy +x. Thus N, is the observed echo 


rate per unit length of echo line. 


* If the height interval from which echoes are obtained is reasonably small compared with the 


mean height h we can treat the minimum detectable line density as constant over the strip and equal 
to its value ap on the mean echo line. 
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If O(«,) can be represented by equation (8) we get the simple result: 


N,= Al o( - ) (14) 


siny \cosx 
where 


fae ge, ese . (27) * 86-1 (e—1) 
I= (2) | 2(1 —u?) tdu= (=2) T'(s—2) (1 


and is plotted in Fig. 7. 
H1 is simply the area of the normalized distribution in height of the reflecting 
points in the strip of the echo plane and J is identical with the integral 


00 a» X(s—1) 

3¢ ' 

I, =| | si; dx 
= 2: 77s 


which appears in earlier work (1, 2, 5,6). HJ/cosy may be regarded as an 
average effective length of the trails of radio-meteors. 











: 27\8-1 71, 7 
Fic. 7.—I= { — | u°-? (1 —u?)*-! du. 
0 


4 


The zenith distance x in the above is the value at the point P on the echo 
line and it differs from the value yo at the observer. To a close approximation 
we may substitute (see Appendix) : 

cosx~Cosyxp, sin? ysin? y,+2h/R, 
where R, is the Earth’s radius, whence we obtain 
N, = Af, (xo, s)O(ap) (16) 


where 
fi(xo s) =I (cos xo)*— (sin? y+ 2h, R,)" ., (17) 


f,(xo,5) is given in Fig. 8 as a function of yx, for a range of values of s (with 
= 100 km). 
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Using the method of Section 2, the sensitivity factor p may be determined at 
all points along the echo line, hence the echo rate N over part or the whole of the 
echo line can be obtained, by numerical integration, from 


N= Hf xoy5)0(¢0) | p*1 dl (18) 


The echo rate between specific range limits can be found by evaluating the integral 
in (18) along the appropriate intercept of the echo line. 


20 











0 





Fic. 8.—The rate factor £,(xo, s) as a function of x» for various s (h=100 km). 


The value of s can be found from the distribution in amplitude of the shower 
meteor echoes or from the dependence of echo rate on equipment sensitivity 
(5,6, 15). By substituting it into (18) we obtain the factor relating N to 0 
and hence the incident flux may be deduced from the observed echo rate. 

5. The echo rate per unit range interval.—In Fig. 1, @ is the angle between OP 
and the normal to the echo line at P. ¢, ¢, are the elevation of OP above a 
horizontal plane at P and O respectively. 6 is the angle QOP where OQ is the 
intersection between the echo plane and the plane through O, the zenith and 
radiant point. 

We have sin@=dR/di where dR is the range interval corresponding to the 
element di of the echo line. Thus if N,dR is the rate of echoes with ranges 
between R and R+dR we obtain, using (14), 


N= Ht o(22-). os) 


~ sin @sin x \cosx 





To a sufficient order of approximation (see Appendix), 


sin 6 sin x = +/ (sin* xp — sin* do), COS x COS Xp; 
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hence 
Ny= Hfg(x0,5)®(#p) (20) 

where 
fo(Xo» 8) = J (cos xo)°~* (sin® Xo — sin® Jy)~**; (21) 

N, depends on R through sin ¢y ~ : _ x 

Equation (20) gives the rate per unit range; it can be applied to a finite range 
interval provided the sensitivity factor p does not vary appreciably along the 
intercept of echo line. It clearly fails when ¢, approaches y,, when the echo 
line becomes tangent to the circle of constant range R (i.e. when P and Q in 
Fig. 1 coincide). In this case we must revert to the method of Section 4 and 
evaluate the integral of p*-! along the echo line (for the result to be valid the 
range interval must be sufficiently large to encompass the range spread due to 
the finite dispersion in height of the reflection points). In this manner the mean 
rate N’ between 400 and 500 km has been evaluated, for echo lines intersecting 
the beam axis of the Sheffield rhombics at range 450km*. The result is 
illustrated in Fig. 9 where the full curves give N/[HO(ap.)]=fe(xo,5); they 
may be compared with f,(y9, 5) evaluated for R = 450 km, given as broken curves. 





20 











——_—$f, (xo, 8) evaluated for the Sheffield rhombic aerials with 400< R<500 km. 
———fa (xo, 5) for R= 450 km ($9= 10°9°). 


6. The determination of shower radiants.—The method of Clegg (7) and its 
modification due to Keay (8) rests essentially on the assumption that the 


* This will be approximately the condition for maximum rate per unit range at R=450km 
(but see Section 6). 
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maximum shower echo rate at a given range occurs when the echo line intersects 
the beam axis at that range. Referring to Fig. 5 we see that the radiant 
declination and right ascension can be estimated from the times of maximum 
rate on two aerials at different azimuths. It is clear from the above, however, 
that the zenith distance will be changing as the echo line moves through the beam 
causing a displacement in the time of maximum rate which will occur earlier or 
later depending upon whether y, is an increasing or decreasing function of T. 
The magnitude of the displacement will depend, inter alia, on the azimuthal 
width of the aerial beams (being greater, the wider the beam) and on the radiant 
declination. It is clear that the present theory enables these time displacements 
to be estimated and appropriate corrections to be made. 

7. An alternative approach to the problem.—If we consider again a narrow 
strip in the echo plane, normal to the echo line at P, we see that a meteor trail 
with maximum line density x can be detected over a height range Hn (u,/ug) 
where u,, U, (u,> uy) are the positive solutions of the cubic equation (10) with 
%=%p. The quantity In (u,/u.) is thus a function of «m (for given ap) and may 
be regarded as a weighting factor which when multiplied by the differential 
distribution function ®(am/cosy) gives the distribution of maximum line 
densities of the meteor trails which are actually observed (see Fig. 10). 


mame UR (u, u,) 





A ht L i i 
04 05 06_, 07 08 
<_< (~m /%p) 





Fic. 10 


Thus the number of echoes observed per unit length of echo line at P from 
meteors with line densities between %m and am+dam is n(%m)dam where 


alee ey eo(S). (22) 


sinxcosx Uz, \cosx 


For the special case when ® may be represented by the power law we get 
n(am)o¢ 2~*In (u/s) (23) 


where =%m/ap. When z>1, u,/u.~272/4. Fig. 11 illustrates how n(#m) 
from (23) varies with am. For s=2 the most frequently observed meteor has 
maximum line density about 1-3 times the threshold value, «». 
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The total observed rate per unit length of echo line, for any distribution ® 
is obtained by integrating (22) giving 
, H > u skp 
Biestie “In “1 at) de (24) 
sin x COS x u, \cosyx 
Equation (24) enables the observed rate to be predicted if ® is known but 
unfortunately it is of little value in deducing ® from the observations. It is 
nevertheless instructive in showing that only for the power law does N, have 
the same functional dependence as © on the limiting line density, «,; in this 
case (24) becomes 


N,= (s—1)H 4 o( Xp ) [“s-ein (u,/u,) dz. 
rae | 


sin x cos x 








arbitrary units 


4 





rt 
6 
_e Xm/~p 





Fic. 11.— Differential distribution in maximum electron line density of observed meteors for an 
incident distribution satisfying an inverse power law with exponent s. 


8. Discussion.—A number of approximations have been made in order to 
secure the simplest procedure for deducing the incident flux of meteoroids 
from the observed radio-echo rates. ‘To whatever extent this may limit the 
absolute accuracy of the result, the above treatment should prove of considerable 
value for comparison of the incident fluxes and magnitude distributions in the 
various meteor showers. 

Errors will clearly result if the ablation process deviates from the theoretical 
model; for instance if the phenomenon of fragmentation (14) is important 
for the faint radio-meteors the factor HJ in the echo rate equations will be an 
overestimate. Rapid diffusion of meteor trails will limit the height to which they 
can be observed, particularly on short wave-lengths, while on longer wave-lengths 
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(A> 15m) ionospheric D-region absorption and Faraday rotation may at times have 
a significant effect on the rate (Baldwin and Kaiser (unpublished) have observed 
a marked decrease in the echo rate on 17 Mc/s associated with a short-wave radio 
fade-out). The effects of rapid diffusion will be most marked for fast meteors 
(which tend to ionize at greater heights) and for showers with values of s markedly 
less than 2, when the theory predicts a substantial proportion of echoes from 
meteors considerably brighter than the threshold for detection (and which may 
therefore be seen to greater heights). 

Deviation of the distribution of maximum line densities from the simple 
power law with constant exponent s will lead to error since it was shown in 
Section 7 that only in this case does the echo rate have the same functional 
dependence as © on a». Observations of sporadic meteors made at Sheffield 
(unpublished) give s near to 2-0 down to 11th radio magnitude, with no significant 
seasonal or diurnal variation. However, considerable deviations from a power 
law are found in the case of a number of meteor showers. In these cases the 
exponent deduced from the variation of echo rate with equipment sensitivity 
(or from the echo amplitude distribution) will be close to the true value only 
if s varies sufficiently slowly with «,, this condition being more stringent the 
smaller the value of s. In fact, provided s is not too small, most observed meteors 
will have magnitudes near to the limit of detection in which case the threshold 
value of s may be used in the echo rate equations without serious error. 

It has been assumed that the echoes ail originate in a narrow band about a 
mean height h (which will tend to be higher the larger the meteor velocity and 
the more sensitive the apparatus) and that the limiting sensitivity is constant 
across this band. In fact the minimum detectable line density will tend to 
increase with height across the echo band (due to increased range) and this will 
again cause the present theory slightly to overestimate the echo rate. The 
representation of the echo band by a mean echo line on a constant height surface 
introduces an additional limitation since the echoes corresponding to the point P 
on this line will actually be spread in range over an interval AR where 


AR ~ HI (sin¢d)—! = HI (sin? dg + 2h/R,)-*. 


This will not be serious provided (a) the sensitivity factor p does not vary 
substantially over the range interval AR and (bd) that if x, approaches ¢, the range 
interval within which the echoes are counted is at least a few times AR. 

It has been further supposed that the shower has a well defined point radiant. 
A spread in the radiant will not be serious provided it is reasonably small compared 
with the aerial beam width, except when yp approaches ¢, when the range interval 
needs to be large enough to encompass any additional range spread due to the 
finite radiant area. 

A final difficulty can arise in subtracting the sporadic rate from the total in 
order to obtain the shower rate. Since s is generally smaller for shower than for 
sporadic meteors the ratio (shower rate: sporadic rate) decreases with 
increasing equipment sensitivity (increasing total rate) until eventually the 
shower may be scarcely observed against the sporadic background. In these 
cases it may be difficult to estimate shower rates unless some method for 
selecting shower meteors (e.g. by velocity) is used to distinguish them from 
the background. 
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APPENDIX 


(i) Transformation of the radiant coordinates. 
Let 
1=latitude of observer, 
8=radiant declination, 
T= radiant hour angle, 
gr = radiant azimuth (E of N), 
Xo = radiant zenith distance. 
A convenient form for the solution to the astronomical triangle is 
f(xo)—sin Xo(1 — cos Yr) = f(5 —/) (26) 
where 
f(x)=sinx+tan/cos x. 

Given the latitude /, f(x) may be presented in tabular or graphical form and 
used to evaluate 6 for any given yo, gr. The hour angle can then be obtained 
from either of the relations 

SiN Xo SIN yr 


sin T= . 
coso 


(27) 


COS x) — sin/sind 


cos 7'= - 
cos/cosé 





(28) 


(ii) Equations for the echo line 
If do, g are the elevation and azimuth from O of the point P on the echo line 
(Fig. 1) then it can be shown that 
tan dy= — tan x, cos (¢r— @). (29) 


The slant range R and the curved distance Rs (see Fig. 12) are related to dy 
through 


R~R, (,/sint y+ . —sin $0) (30) 
E 


and 
Rs~= Roos d, (31) 


where R,, is the Earth’s radius. 

Equations (29), (30), (31) give the coordinates (Rs, gy) of the point P on 
the echo line for any given zenith distance x». They have been used for the 
construction of Fig. 3. 


(iii) Equations for the curves of Fig. 5 

We wish to determine the radiant hour angle when the echo line corresponding 
to a radiant at declination 6 intersects a point on the echo surface at azimuth 
and slant range R. R is related to ¢y by (30) and @ is obtained as a function 


of ¢) from (29). 4 is then obtained as a function of JT through the procedure 
in (i) (note that yo >¢p). 
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(iv) Relations between $, do, x, Xo, 9 (see Sections 4, 5) 


Consider the triangle OPC (Fig. 12) where C is the centre of the Earth. 
It will be recalled that ¢ is the elevation of OP above the horizontal at P, thus 
the angle OPC is equal to 7/2—¢ and we obtain 
en ee 
from which it follows that 
2h 


sin? ¢ = sin* dy + 7 (33) 





Fic. 12 


In the right-handed coordinate system Oxz of Fig. 12, let the radiant have 
direction cosines (/, m,n) where n=cos x9. The direction cosines of OP and CP 
are (cos¢p, 0, Sindy) and (sinc, 0, cosc) respectively. Now x is the radiant 
zenith distance at P, hence 

cos x =/sinc¢ + cos x9 Cos ¢ 
and since the radiant direction is normal to OP, 
0=/cos by + COS Xp SIN do. 
Eliminating / from the above two equations, and putting 6=¢)+¢, we get 
cos cos h 
Ie wn A ns +—I 
cosy  cos¢d Ry 





(34) 
: ' 2h 
and sin? y= sin? yy + —. 
Ry 
Referring to Fig. 1 we see that sinycos#=sing, whence, using (34), we 
obtain 


sin x sin 6 = (: + Rk) / (sin? x9 — sin? do). (35) 
E 


The University, 
Sheffield: 
1960 April. 
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OSCILLATOR STRENGTHS OF NEUTRAL ATOMS 
OF THE IRON-GROUP 


C. W. Allen 
(Received 1960 May 11) 


Summary 


Measurements from various sources of the relative oscillator strengths 
of multiplets in Sc, Ti, V, Cr, Mn, Fe, Co, Ni have all been reduced to a 
consistent absolute scale. This makes it possible to compare experimental 
accuracy, make certain corrections, and tabulate mean measurements. The 
measurements are then compared with calculations based mainly on an 
elaboration of the f-sum rule. The character of the results differs for three 
sets of multiplets as follows :— 

(a) For the transition types A, E, F, G, H, which involve non-equivalent 
electrons only, the log (observed/calculated) intensity is almost independent 
of excitation potential or wave-length. ‘The scatter is probably due to 
unsuitability of the designations or to the interaction of terms. When 
absolute f-value measurements (for Cr, Mn, Fe, Co, Ni only) are used to 
convert these multiplets to a measured absolute scale, all atoms are within a 
factor 3 of the calculated scale. The high (obs./calc.) value in Mn is 
discussed. There is only a small systematic difference between observation 
and calculation. 


(5) For the transition types B, C, which involve equivalent electrons, 
the scatter of the results is much larger. Miultiplets with an upper 
excitation potential less than 4 eV have f-values below the calculations, 
while above 4 eV values agree with the calculated absolute scale. 

(c) Measured values of multiplets violating LS-coupling vary erratically 
between zero and an empirically calculated value. 

Recommended corrections for putting the mean values onto the best 
absolute scale are given. 





1. Introduction.—The atoms that give rise to the greatest number of absorption 
lines in stellar spectra are mainly in the first long period of the elements. Of these, 
Sc, Ti, V, Cr, Mn, Fe, Co, and Ni have rather similar spectroscopic properties 
typified by their most abundant member Fe. They may be called the iron-group 
atoms. Since these elements have such rich spectra the oscillator strengths of 
their lines are in strong demand for astrophysical analysis. Unfortunately 
oscillator strength measurements of these atoms on an absolute scale are 
attended by considerable difficulty and possible error. Moreover the atoms are 
too complex to allow reliable calculations. Both measurements and calculations 
are being used and it is not certain which gives the best absolute values. 

The object of the present work is to compare the experimental and calculated 
oscillator strengths in the iron-group atoms. We thus obtain an indication of what 
regularities and general trends occur and what values are reliable. 
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Although measurements and calculations may be applied to individual 
spectrum lines the number of lines concerned is too great for convenient 
tabulation and analysis. The situation may be relieved by using data for complete 
multiplets only. The values of the weighted oscillator strengths gf, where g is 
the statistical weight and f the oscillator strength, may readily be added within 
each multiplet to give the weighted oscillator strength of the multiplet g;/f. 
The subscript t implies that the values refer to complete terms and multiplets, 
not levels and lines. The symbol gf must refer to the multiplet as a whole 
since f is not the same for each line. When summing gf to give g¢f an allowance 
can be made for missing or unreliable lines by using the rules for intensities 
of lines within multiplets (Condon and Shortley 1935, p. 241; Allen 1955, p. 55). 
Usually the best measured or calculated lines are the strong ones and the 
contribution from the rest of the multiplet does not give rise to much ambiguity. 
If the intensity rules do not apply, as for multiplets forbidden in LS-coupling, 
the eye estimates of intensity in the Revised Multiplet Table (Moore 1945) 
usually lead to a suitable estimate of the conversion from gf to g:f. In the 
present paper all values are expressed as g¢f or log g¢f. 

2. Measurements.—Some relative oscillator strength measurements have been 
made for all of the eight iron-group atoms. ‘There are satisfactory absolute 
measurements in Cr, Mn, and Fe, and also some absolute values have been 
obtained for Co and Ni by quantitative comparisons made in arc spectra of 
diluted copper alloys. Only relative values are available for Sc, Ti, and V. 

In order to compare and average the experimental relative gf values it is 
important that they should be systematically reduced to the same scale. The 
ideal would be to convert them all to an experimental absolute scale, but for some 
atoms no absolute values exist and for the others the absolute values are likely 
to change as better measurements become available. In order to achieve some 
stability and consistency 1 have converted all the measured relative values to a 
calculated absolute scale. Simple methods of obtaining calculated values are 
given in Section 3. 

First the experimental values have been converted to log gf in the units of 
the published experimental work. Where the results are quoted as intensity / 
the value used is A*J with A in microns, or where transition probability 4 is 
quoted we use A2gA. Sometimes a Boltzmann factor or a correction to a Boltzmann 
factor has been applied. A reduction factor has then been found by plotting 
experimental against calculated values. In doing so, greater weight was given 
to the higher values of gf and to the types of transition found to give more regular 
results (see Section 4). The reduction factors and any Boltzmann corrections 
or factors are given in Tables | and II. 

Both measured and calculated values of g;f are set out in Table 1. The 
columns are: the Multiplet Number from the ‘‘ Revised Multiplet Table ”’ 
(Moore 1945); the designations of the lower and upper terms from ‘* Atomic 
Energy Levels’’ (Moore 1949, 1952); the type of transition using the notation 
of Section 3; the lower and upper excitation potentials; the wave-lengths A 
of the leading lines; relative measurements of logg;f reduced to a calculated 
absolute scale (the asterisk * implies a temperature adjustment has been made) ; 
mean values of logg;f on the same scale (the experimental results have been 
used almost entirely but a little weight has been given to calculated values in a 
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few cases where this seemed appropriate); calculated values of logg:f from the 
f-sum rule (Section 3); and calculated values from the Bates-Damgaard (1949) 
procedure. Except for a possible absolute scale adjustment the ‘‘mean”’ 
column is intended as the best available estimate of true logg¢f. 

Sources of experimental values for various atoms and the author abbreviations 


are as follows :— 


KK King and King (1938) Ti, Fe 

King King (1947, 1948) 

HK Hill and King (1951) Cr 

KPDO King, Parnes, Davis and Olsen (1955) Co 

OPP Ostrovsky, Parchevsky and Penkin (1956) Ti, Mn 

OP Ostrovsky and Penkin (1957 a, 1957 b) Cr, Sc 

Ost Ostrovsky (1958) V, Co 

AA Allen and Asaad (1957) Cr, Mn, Fe, Co, Ni 
Mit Mitrofanova (1952, 1954, 1955, 1958) Fe, Ce, Ti, Se 
vSS van Stekelenburg and Smit (1948) Ti 

AHB Aarts, Harting and Bakker (1954) Fe 

Car Carter (1949) Fe 

vD van Driel (1935). Data from King (1948) Ni 

OB Ornstein and Bouma (1930) Co, Ni 

HD Heid and Dieke (1954) Ni 


The more accurate and reliable values of loggtf are quoted to two places of 
decimals. For the ‘‘mean’’ column, reliability is determined by interagreement 
of experimental results or agreement with calculations. 

All measurements of oscillator strengths require a Boltzmann factor to allow 
for the declining population of higher levels. The factor may be very large and 
may cause serious error if the excitation temperature 7’ is not well determined. 
The factor is usually applied in the form 


log gf (final) = log gt f (as measured) + E(5040/T) 


where E is excitation potential in volts. If the experimentally determined 
(5040/7) is unsuitable the final logg;f results require a correction E A(5040/T). 
Corrections have been made for some atoms (labelled *) in Table I. In order 
to keep the temperature correction and scale correction as independent as possible 
the temperature correction has been pivoted on a representative excitation 
potential for which we adopt 1 volt for E, (lower E.P.) or 4 volts for E, (upper 
E.P.). Thus the corrections to log g¢f are 


for absorption lines A(5040/T)(E, — 1°0), 
for emission lines A(5040/T)(E,— 4:0). 
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The corrections and factors in Table II have been applied in Table I and are 
discussed below or in Section 4. 


Tase II 
Boltzmann corrections and factors 





Cr HK* 
Cr AA* 
Cr Mit* 


A(5040/T) 





0°30 
—0'07 
—0'50 


A(5040/T) 





—0°05 
—0°20 


Co OB* 
Ni OB* 
Ni HD* 


(5040/7) 





1°00 
0°95 
1°00 
































The differential temperature corrections A(5040/T) and the experimental 
accuracy may both be studied by plotting observer comparisons against excitation 
potential. In Fig. 1 we give a plot of log(g¢f 1st observer/g¢f 2nd observer) 
against E.P. for those cases in Table I where a number of multiplets have been 
observed by several authors. An upward slope means that the 1st observer 
has used a temperature that is low on the standard represented by the 2nd observer. 
We note for example for Ti that the OPP temperature is low by comparison 
with KK whereas vSS and KK agree. More details will be considered in 
Section 4 when the results are compared with calculations. 
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Fic. 1.—Comparison of observers ploited against lower excitation potential. 
The logarithmic scale is indicated by 1 dex (=log,, 10). 


Several of the experimental excitation temperatures have been derived by 
comparison with the work of King and his collaborators and therefore agreement 
of slope in Fig. 1 does not necessarily mean independent temperature confirmation. 
The measurements by Mitrofanova appear to be seriously influenced by self- 
absorption. Without allowing for this her excitation temperatures are derived 
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by comparison with King. This procedure apparently leads to large errors in 
the derived values of g:f and a further empirical temperature correction is 
sometimes required to make her results fit other work and become usable. Some 
of the Allen and Asaad results have been subject to some revision which has 
resulted in minor changes of temperature. 

The scatter of the points in Fig. 1 gives an indication of the accuracy of the 
experimental work, and leads to the approximate errors in logg:f given in 
Table III. 


Tas_e III 
Approximate standard errors in log gtf 








KK, King, HK, KPDO +009 OB +0°25 
OPP, OP, Ost +o°12 vD +0°18 
AA +0°15 AHB +0°13 
Mit +0°5 Car +0°18 
vSS +0°08 HD = +018 




















3. Calculations.—Two approaches to the calculation of oscillator strengths 
that have some prospect of general applicability in the iron-group atoms are 
(a) the f-sum rule, and (5) the wave-function method systematized for a coulomb 
approximation by the Bates-Damgaard (1949) procedure. Since we cannot 
expect results to be more accurate than about +o-1 dex (i.e. +01 in the 
10-based logarithm) we are entitled to introduce various simplifying approxima- 
tions. 


f-sum rule 
The f-sum rule (Unsdéld 1955, p. 350), written in terms of gyf becomes 


Labi f+ Dei ft Dsrfl=re (1) 


where t represents the atomic term for which the summation is being 
made and g; is the total weight of that term; u signifies transitions to upper 
terms, c to the continuum and | to lower terms; r is the number of electrons 
involved in the transition. The summation }ig:f is an emission oscillator 
strength and therefore negative. The other summations include transitions 
to high series members and to the continuum. However, the spectrum lines 
that concern us in the iron-group atoms are mainly lowest members of their 
particular series and as such they contain most of the total oscillator strength. 
We suppose the oscillator strength of the lowest member to be 5b times the total 
of upward transitions (including the continuum), where b<1, then 


Dabif =brg,— Digf (2) 


where the summation on the left is now for the lowest member in each series and 
is usually over one or two arrays (hence subscript a). We adopt b=o-7 on the 
argument that in simpler atoms where 5 can be determined it is between 0-4 and 
10. Later we compare this adopted value of 6 with both Bates~Damgaard 
calculations and experiment. The emission component —)igtf is often zero 
for the arrays with which we are concerned; otherwise it is positive. Apart 
from sign the numerical values of g¢f are the same in emission and absorption. 


24* 
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The expression (2) enables us to compute a value of gf (though not 
necessarily a sound one) for nearly all multiplets that are amenable to any sort 
of calculation. The formula is of very general use when applied to iron-group 
atoms. It is convenient to distinguish the atoms by z as follows: 


Sc Ti V Cr Mn Fe Co Ni 
n= 2 3 4 5 6 7 8 9 


Results of calculations are often expressed in relation to the weights of various 
terms, g,=(2S,+1) (2£,+1) where 2S,+1 is the multiplicity and L, is the total 
orbital quantum number for the term t. 

We give below the notation and the method of calculating gf for the various 
transitions : 

A=d"s—d" p (examples in all atoms Sc to Ni); r=1 since only 1 non-equivalent 
electron moves; }i=0; 6b=o0-7. All transitions are from lower unit 
terms to upper triads with three times the total weight of the lower term. 
&tf is proportional to gut the weights of the upper terms. From (2) we have 


&if=0°78ut/3 =0'2338ut- 

B=d"-1s*-d"-1sp (examples in all atoms Sc to Ni); r=2; Si=0; b=0'7. 
Each term t combines with two upper triads. For one of these the parent 
multiplicity is (2S,+1)+1 and for the other (2S;+1)—1. Assuming the 
gtf to be proportional to the parent multiplicity (see Rohrlich 1948) we 
obtain 

&t f=0°233gut(2S¢+1 + 1)/(2S¢+1) 
where the + and — refer to the higher and lower parent multiplicities. 
The number following B in Table I is (2S;+1+1). 

C=d"s—d"-' sp (examples in Ti, V, Cr, Fe, Co, Ni); r=n; Si=0; b=0-7. 
To apply the f-sum rule for the moving d electron one must add the 
transitions to both d"-'sp and d"-1sf which are in the ratio 1:3} (see 
Allen 1955, p. 60). The total gtf in the d"s-d"-!sp array will then be 
0-7(2/9)ng:. Comparing this with c/(.#), the sum of the angular matrices 
as published by Goldberg (1935) with the correcting factors (Goldberg 
1936, Allen 1955, p. 60), we obtain the generalization 

&tf=0-7e% (M)/27 =0°02607(.#). 
When the c¥(.#) for the individual multiplet is not given, but only the 
sum of x multiplets, it has been assumed that the multiplets actually 
concerned are stronger than the average. To cater for this effect I have 
divided the total by 1/x; not by x. Large errors in the calculations are 
inevitable. 

D=d"*!-d"p (one example in Ni); r=10; }i=0; b=0°7; (2/9) of total gf 
as in C, 

E=d"~-!sp-d"-1 ss (examples in V, Cr, Mn, Fe, Co, Ni), 

F =d"~-!sp-d"— sd (examples in Ti, V, Mn, Fe, Co, Ni), 

G=d" p-d"s (example in Cr), 

H=d"p-d" d (examples in Cr, Ni); r=1; — Dig, f~0-233g,; b=0'7. The ratio 
of the p-s to p—d transitions is 1 to 10. Hence the /-sum rule gives 

EandG = Dg,f=(1/11)(0-7g,+ 02338), 
FandH 5g,f=(10/11)(0-7g,+0°233g;)- 
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Comparing these values with c/(.#) as obtainable from Goldberg’s tables 
(1935, 1936) we obtain for E, F, G, H 
&,f =0°08507 (M). 

J =d"—' sp-d"— p? (one example in Ti); r=2; Sig,f~0-233g,; b=0-7. The 
f-sum rule (2) gives >, 9,f=0'7 x 2g,+0°233g, Comparing such values 
with the tabulations of c/(.#) we obtain 

&,f=0°817e/ (4). 

h The symbol h is added to transitions which are not the lowest members of 
the series and to which therefore the above rules do not apply. B.D. 
calculations are used in such cases. 

In addition to the transitions permitted in LS-coupling there are many strong 
transitions that violate the LS rules. The following notation is used: 

ic =intercombination line, 

ip =interparent transition, 

ml= multiple change in L, 

dj =double jump. The main such transitions are of the forms d"~!s*-d"p, 
d"—1 sp-d" d, d" p-d"—" ss, 

In such cases, the f-sum rule can give an approximate maximum value. By 

comparison with permitted multiplets it would appear reasonable to adopt for 
our calculations g,f<0-233g, where we select t to be the lower term of the 
transition. 


Calculated values of g,f in accordance with these f-sum applications are 
given in the f-sum columns of Table I. 


Bates—Damgaard calculations 


The Bates-Damgaard (1949) method of calculating the radial integral 
designated o? may be used for many of the multiplets permitted in L.S-coupling. 
The multiplet strength is then given by o?c%(.#) and the weighted oscillator 
strength by 

8, f =0°0304070/ (.M)/A 
with o? and c%(.#) in atomic units, and A in microns. 

The B.D. results quoted in Table I have been taken from any available 
calculations (many from Asaad 1956) and may not be entirely consistent. 
They are not required for the main part of the present analysis. 


Comparison of f-sum and B.D. calculations 


The f-sum and B.D. calculations are compared logarithmically in Fig. 2. 
For the A transitions the two calculations would have been almost identical 
if 6 had been taken as 0-82 instead of 0-70. For the B transitions there is 
considerable spread between the two calculations; the choice of b=0-76 would 
have produced slightly better mean agreement. For the E, F, G, H transitions 
the two calculations are in good proportion and absolute agreement would 
have been obtained by selecting b=o-50. Our choice of b=0-70 gives fair 
overall agreement between f-sum and B.D. calculations. The disagreement in 


calculated values is small by comparison with the disagreement between 
observation and calculation. 
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Fic. 2.—Comparison of log g,f calculated from the f-sum rule and from the Bates-Damgaard 
(1949) procedure. Transition types are A; B; and E, F, G, H. 


It should be noticed that the agreement between f-sum and B.D. calculations 
is valid over a considerable wave-length range. The B.D. procedure has been 
recomputed for each line or multiplet and should not produce an error that 
changes systematically with wave-length. Hence when a wide wave-length 
range is to be considered, sum rules should be applied in terms of f-values, as 
has been done here—not in terms of line strengths which are proportional to gfA. 

To represent the calculated values I have chosen the f-sum procedure because 
it is very general, quite consistent, and easy to compute. It agrees slightly 
better with unadjusted experimental values than the B.D. procedure. 


4. Comparison of observations and calculations.—In Fig. 3 the observed values 
taken from the ‘‘mean’’ column of Table I (and therefore adjusted to a 
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Fic. 3.—Comparison of measured and calculated values of the weighted oscillator strength. 
Arrowed signs (<, <-, etc.) are violations of LS-coupling. 
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calculated scale) are plotted logarithmically against the adopted (f-sum) 
calculations. The transition types are sub-divided by symbols. Arrowed signs 
(e.g. <, <, etc.) indicate that LS-rules are violated and that calculations can 
give an indication of only the maximum value. Such multiplets are included 
in Figs. 3, 4, 5 and 6 only if the log (observed/calculated) value is greater than 
—1'0. 

In considering Fig. 3 and later diagrams, it is well to remember that even 
for genuine regularities several exceptions may be expected. These may be due 
to interactions between terms and to unsuitable designations. Such intensity 
effects may be large and one could not expect the consequent scatter to follow 
random laws. We will therefore look for regularities that relate to the majority, 
but not all, of the multiplets. The multiplets that are strong both by observation 
and calculation are likely to be well designated and unaffected by interaction. 
These, which will appear on the top-right of each plot in Fig. 3, should show the 
best agreement. 


The agreement shown in Fig. 3 varies from atom to atom but we can draw 
the following generalizations : 
(i) The transitions A, E, F, G, H, involving the movement of a non- 
equivalent electron, show observations and calculations in fair agreement. 
The absolute scale of the observations has been adjusted to make this 
agreement as close as possible. 
(ii) The transitions B, C, involving the movement of an equivalent electron, 


show much larger scatter and the observed values are systematically smaller than 
the calculations. 


(iii) The most intense of the observed multiplets that violate L.S-coupling 
are in near agreement with calculated maximum value. 

In Fig. 4 the log (observed/calculated) values for each atom are plotted 
against the lower excitation potentials, E,. We see that for most atoms the A, E, 
F, G, H transitions (non-equivalent electrons) do not change systematically 
with excitation potential (E.P.), while for the B, C transitions (equivalent 
electrons) the low voltage multiplets have low observed intensities. To show 
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Fic. 4.—Relation between the (observed|calculated) gf ratio and E, the lower excitation 
potential. Broken lines show possible trends. Arrowed signs (A, +, etc.) are violations of 
LS-coupling. 
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the trends more clearly, Fig. 5 collects for all atoms (i) the A, E, F, G, H 
transitions, (ii) the B, C transitions, and (iii) the more intense of those multiplets 
that violate LS-coupling. In Fig. 5 (i) there appears to be no systematic departure 
from the zero ordinate (i.e. obs.—calc. agreement). However, in interpreting this 
diagram we should remember that the absolute values have been chosen to render 
good absolute agreement and this may have some effect on the apparent change 
with E.P. Fig. 5 (ii) shows that the characteristics of the B and C transitions 
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Fic. 5.—General relation between the (observed|calculated) g,f ratio and E, the lower 
excitation potential. Broken lines show possible trends (see text). Symbols as in Fig. 4. 


are somewhat similar. They both have systematically low observed values at 
low E.P. but tend towards agreement with calculations for E,>2eV. We see 
later, Fig. 7, that the scatter of the values is somewhat reduced if we use the 
upper E.P., E;, instead of the lower. From Fig. 5 (iii) there appears to be an 


upward trend of the highest observed—calculated values for multiplets forbidden 
in LS-coupling. 

In Fig. 6 the log (observed/calculated) values for all atoms are plotted 
against wavelength. There appears to be no significant trend for A, E, F, G, H 
transitions while the trends for B and C are somewhat similar. Any apparent 
trend for the multiplets violating L.S-coupling is probably not significant. 
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Fic. 6.—Relation between the (observed/calculated) g,f ratio and wave-length. Broken lines 
show possible trends. Symbols as in Fig. 4. 


We are now in a position (ref. Figs. 1 and 4) to compare the measured and 
calculated values for the individual atoms and to investigate what effects might 
have been caused by the temperature adjustments that have been made to the 
experimental results. 
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The A transitions agree while the B disagreement is large and typical. 
No temperature adjustment has been made to the OP observations and 
there is no evidence that any is necessary. If x*P° is designated (from 
Sc11) as configuration 4s(4S)4p the calculated logg:f becomes ov1g in 
agreement with observation. 

The A, F and J transitions are mostly in agreement and the B disagreement 
is typical. Some C measurements are high. The pattern of agreement 
is not greatly changed if the Rohrlich (1948) classification is introduced 
instead of that tabulated by Moore (1949). The upward trend of the 
A transitions in Fig. 4 suggests that a higher temperature would improve 
the agreement. However the results are dominated by the KK observa- 
tions and Fig. 1 shows that the OPP temperature standard is lower still. 
There is therefore no experimental justification for making an upward 
temperature adjustment. 

The agreement is good for the strong multiplets of the A, E, F types but 
the weaker A measurements are systematically small. Fig. 4 shows that 
no temperature change would improve the agreement and therefore this 
atom confirms the King temperature scale. Fig. 4 shows no upward 
trend in the B multiplets which is unusual. 

The general disagreement between measurement and calculation in Cr 
is too large for any precise conclusions to be drawn. Fig. 4 shows an 
upward slope if the highest values are used but this may be associated 
with the general scatter. In any case Fig. 4 implies that the upward 
temperature correction given to the HK measurements was necessary, 
and furthermore Fig. 1 shows that the OP temperature scale is higher 
than HK®* (even after correction). There is other evidence (Allen and 
Asaad 1957, pp. 39, 41) that the HK temperature used for Cr was 
abnormally low. 

Agreement for the A, E, F type transitions is surprisingly good for this 
atom. A zero E, Btransition also agrees whichis unusual. No temperature 
adjustment has been made to the AA scale, which is effectively the King 
scale, and none appears necessary. 

The A, E, F transitions are in fair agreement and the C transitions rather 
fainter than in some other atoms. The experimental temperature agree- 
ment shown in Fig. 1 is satisfactory but the temperatures used are not 
independent. The trend in Fig. 4 shows that a higher temperature 
(about —o-13 in 5040/T fitting the broken line) would lead to better 
agreement with calculations. 

The general agreement is rather poor for the A, E, F transitions. The 
other transitions are low and typical. The character of Fig. 4 is rather 
dominated by the selection of 5040/7'= 1-0 for OB* but Fig. 1 shows this 
selection to be in reasonable agreement with KPDO. ‘Therefore the 
absence of any trend in Fig. 4 is probably correct. 

The right-hand side of Fig. 4 has been forced to agree with the left by 
selecting the temperature of OB* and the reduction factor for HD*. 
There is therefore no unbiased experimental evidence for the good general 
agreement that appears in Fig. 3 for the A, E, F, H transitions. More 
experimental work is needed to connect the observations in the violet 
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with those in the green-red. It is unusual for the low-E, C transitions 
to agree, but the low B values are typical. The E, range in Fig. 1 is too 
small for any effective temperature comparison. 


5. Absolute experimental values.—The absolute oscillator strength measure- 
ments that are thought to have useful accuracy for iron-group atoms are from the 
authors (with abbreviations) listed below: 


KW Kopfermann and Wessel (1951) Fe. 
HL Huldt and Lagerqvist (1952) 
Est Estabrook (1952) Cr. 
DRK Davis, Routly and King (1954) Cr. 
King King (1942, 1956). King’s (1956) letter 

amends DRK and King values Cr, Fe, Mn. 
AA Allen and Asaad (1957), giving absolute com- 

parisons but dependent on other absolute 

values Cr, Mn, Fe, Co, Ni. 
Z Ziock (1957) Fe. 
OP Ostrovsky and Penkin (1957 a) Cr, Mn. 
BDKR Bell, Davis, King and Routly (1958, 1959) Fe, Mn. 


We summarize in Table IV the latest values published or available and the 
mean experimental values adopted for the present comparisons. 


TaBLe IV 


Absolute f and log gt f values 





Cr Mult. No. 1 Est. King OP Adopted 
£(4254) 0086 0°055 Or15 
f(4275) 0-069 0°043 ovr! 
f(4290) 0048 0°032 0082 
log gtf +0°15 — 0°04 +0°38 +0°12 





Mn Mult. No. 2 BDKR, King HL OP Adopted 
(4030) 0'061 0062 0°056 
f(4033) 0°034 0°047 0°040 
£(4034) 0°032 0032 0°027 

log gtf —oO'12 —0'07 —0'13 





Fe Mult. No. 5 King BDKR KW Adopted 

£(3720) 0°033, 0°022 0°032 0°043 0036 
log (gt f/ef) 058 0°53 
log et f +0°04 


























In Table V the observed values of Table IV are compared with mean values 
in Table I, which are observed values on a calculated absolute scale. They are 
also compared with the AA values which are all on a consistent scale. 
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TABLE V 
Comparisons of absolute log g,f 





log etf log differences 





Observed Calc. scale AA-calc. 
(Multiplet) Mean of Obs.-AA All 
Table III Table I multiplets 








Cr (1) +o0°12 +0°14 —0'22 
Mn (2) —or'll —0°54 +0°32 
Fe (3) +0'04 —o'10 . +0°10 
Co 


— 0°05 
Ni 


— 0°36 


Mean — 0°04 























The greatest difference quoted in Tables IV or V is 0°43 dex or a factor 
of 2-7. From this general agreement we might infer that the standard error of 
obtaining absolute values by (i) f-sum calculations on several well-designated 
multiplets, (ii) well-executed absolute observations, or (iii) the AA dilute alloy 
procedure, is not likely to be greater than +0°3 dex. We must investigate 
however whether any of the three methods can be substantially more accurate 
than this, and whether the larger differences are due to faults in measurement 
or calculation. 

The systematic difference between observations and _ calculations, 
log (obs./calc.) = +018, is not very large and is due mainly to Mn. We see 
from the last column of Table V that Ni differs in the opposite direction and 
we reach the conclusion that there is no substantial systematic difference between 
absolute measured and calculated values. 

Mn is the only atom for which the log (obs./calc.) is large enough to attempt 
a decision on whether the fault lies in the observations or calculations. For Mn 
three independent absolute observations are in very good agreement and 
furthermore the AA observations give further support to the high measured 
values. On the other hand the relative calculated values agree well with the 
relative observations. We must conclude that, either (a) for all seven of the used 
A, E, F multiplets the f-sum (and also B.D.) calculations chance to be too low 
by about the same amounts, or (5) the relative measurements do not give a 
satisfactory connection between the multiplet 2 and the rest of the spectrum. 
Probably both of these points are involved to some extent. With regard to the 
first we might note in Fig. 4 that for all the relevant Mn multiplets Z, is in the 
2-3eV range. If all of the points were raised by about 0-3 dex (to meet most 
of the discrepancy), they would not be out of keeping with the pattern for Ti, 
Cr, Fe and Co in this E, range. 


6. Discussion.—Rohlich (1948) has made a study of the terms in Ti and finds 
that a redesignation of the odd terms will improve the relative f-values of several 
multiplets. The use of Rohlich’s designations reduces the scatter of the Ti 
points in Figs. 3 and 4 to some extent but does not alter the general pattern. 
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His calculations for inter-parent multiplets are very satisfactory but are not 
included in our analysis because there are no corresponding results for other 
iron-group atoms. His investigation leaves the impression that in Ti there is 
no set of designations that can be regarded as completely correct. 

Without doubt the same designation difficulty applies in other iron-group 
atoms. The experience in Ti makes it reasonable to suppose that interactions 
and designation faults can account for most of the scatter of at least the A, E, F, 
G, H multiplets in Figs. 3 and 4. 

We are left with a need to consider the general trends. 

Fig. 5 (i) gives the impression that there is no general trend of log (obs. /calc.) 
against E, for the A, E, F, G, H multiplets. However when the high absolute 
values of Mn are taken into consideration (see end of Section 5) seven points 
in the range E,=2-3eV may need to be raised by about 0-3 dex. Also the 
level of the cluster of points due to Ni near E,=3-6eV has not been determined 
experimentally. It is therefore still possible that there is a vestige of E.P. 
dependence among the A, E, F, G, H multiplets. It could not be more than 
about +0-2dex/eV and therefore very much less than the effect first exhibited 
by Allen and Asaad (1955). Experiments with higher temperature sources and 
independent temperature standards are perhaps the best means of resolving 
final doubts on this question. 

On the other hand the change of the observed/calculated ratio with excitation 
potential is considerable for the B and C multiplets as can be seen from Fig. 5 (ii). 
Also these two transition types are distinguished by having a variation with 
wavelength, as seen in Fig. 6. It seems unlikely that the E and A variations are 
physically independent and there is no correlation between EF, and A that could 
account for the two trends. However, both could be explained by one cause 





o 





' 
nw 





log (observed /calculated ) values 














£, in electron - volts 


Fic. 7.—Relation between (observed/calculated) values (for B and C transitions) 
and the upper excitation potential. The broken line represents the main trend. 


if the departure from agreement were due entirely to the upper level. This does 
correlate with wavelength and could account for the A variation. Furthermore if the 
ratio obs. /calc. be plotted against E,, the upper E.P., see Fig. 7, the relation becomes 
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clearer and the variation larger. This argument suggests fairly strongly that 
it is the position of the upper level that is important in deciding the intensity ; 
that for E,>4eV the B and C intensities agree absolutely with calculations 
(though with a large scatter), and that for E,<4eV intensities are diminished 
by systematic factors that can be as much as 1°5 dex at E,~2-5eV. Although 
this argument points to the position of the upper level as being responsible 
for the depression of intensities of both the B and C transitions there is no evidence 
of a correlation for individual B and C transitions that have the same upper level. 

Finally we should give an indication of what values might be regarded as the 
best available. Except for a possible scale correction these are given in the 
‘mean’ column of Table I. To obtain a scale correction it seems safest to make 
a compromise between observed and calculated values. From the foregoing 
discussion and Table V,I offer the following corrections as the best available 
at present. When no absolute observations are available there seems no reason 
to make any correction. 


Element Sc Ti Vv Cr Mn Fe Co Ni 
Correction in dex oo )600'0—C8Ns—«é*=8lV'LD +0°25 +0°1I0 00 —o'r 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7: 
1960 May. 
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THE PHOTOELECTRIC LIGHT CURVE OF FU ARAE 


G. G. Cillié and E. M. Lindsay 
(Received 1960 May 26) 


Summary 


Photoelectric light curves have been obtained in blue and yellow light for 
the star FU Arae using the 60-inch Rockefeller reflector of the Boyden 
Observatory. There are two equal fairly flat topped maxima. Primary 
minimum is much deeper than secondary; the latter occurs at o°5 phase. 
The light elements are Min. = JD 2426869:490 + 09-864500957E. 





The variability of FU Arae was discovered by Miss Swope (1) in a search 
for variables in Harvard Milky Way Field 233. It was designated Hv 6768 and 
is the star CPD —58°7173. From measurements on 106 plates taken with the 
10-inch Metcalf telescope Miss Swope concluded that the star was an eclipsing 
binary, its magnitude varying between 11-2 and 11-6. She derived the following 
light elements for the system. 


Min. = JD 2426869-490 + 09-86450E. 

FU Arae was observed on 30 nights during 1954 and on 21 nights in 1956 
using the 60-inch Rockefeller reflector and Linnell-King photoelectric photo- 
meter of the Boyden Observatory. The total number of individual observations 
was 1171 in blue and 874 in yellow light. The blue filter was a Corning 
No. 5551; the yellow a Corning No. 3385. Comparisons with a nearby standard 
star are expressed as magnitude differences Amp and Amy for blue and yellow 
light respectively. 

From primary minima observed on two nights in 1956 the following light 
elements were computed: 


Min. = JD 2426869-490 + 09-864500957E. 


Table 1 gives the differences between the observed and predicted minima 
in heliocentrically corrected time. 


TABLE I 
Primary minima of FU Arae 


Observed Min. = 2426869°490 
Date of run time of prim. +01-864500957E 
min. Wt. Calc. min. O-C 
BOD scesen vere 
1956 August 27/28 = 5§713°33335 3 5713°33479 "00144 


1956 September 3/4 5720°25224 3 §720°25080 +:°00144 


In Tables II and III the normal points of the blue and yellow light curve 
are given. The first and second columns give the phases and magnitude 
differences respectively of the normal points. The third column gives the mean 
error for each normal point and the fourth column its weight, that is, the number 
of observations used in forming it. 
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Taste II 


Normal points of blue light curve of FU Arae 


M.e. of 
normal pt. 
m m 
+0°004 05018 +O°112 
*77° 007 5344 *095 
*520 "003 "5614 +0°065 
"004 *5890 —O°'ols 
+0°066 "003 “6120 "022 
003 6361 033 
048 "002 6638 053 
-066 003 “6869 "073 
"003 "7104 "084 
"079 "003 7384 "089 
"094 "005 +7628 "092 
‘090 "002 7882 "092 
"004 ‘S121 
"092 "003 *8370 "059 
"002 8612 038 
‘061 "004 8881 
“005 “QII3 "049 
"003 *9371 
"003 "9631 
005 *9829 *800 


Phase M.e. of 


We. (period) Amp normal pt. 


Amp 


Taste III 
Normal points of yellow light curve of FU Arae 


Phase M.c.of yw, Phase Me. of 
(period) normal pt. * (period) normal pt. 


m m 

0°0000 . +0°'003 34 0°5383 +0°005 

"0272 . 0°02 32 *5858 

"0527 ° "004 28 "6198 

°0837 ° 003 30 6614 . 003 

‘1131 ° "003 24 6856 

°1372 ‘ "002 29 “7108 

"1614 , 003 28 *7392 

*1940 ° "003 25 *7626 

"2251 “ "004 18 *7879 

+2683 , ‘ 22 *8105 

"2970 “ ° 23 *8360 

3246 , , 25 8613 

*3562 ° . 28 +8887 

"3927 ° . 24 "9145 

"4248 , ° 26 "9425 

"4583 ‘ ‘ 38 "9728 

*5018 ° ° 40 


Fig. 1 gives the blue light curve and Fig. 2 the yellow light curve of the star. 
The maxima are of equal height and to a small degree are flat topped. They 
are more so in blue than in yellow light. It would appear that the maximum 
preceding secondary minimum is flatter than that following secondary minimum, 
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particularly in blue light. A weakness in the observational material between 
phases 0-6 and o-7 would result in such an apparent effect but there is no actual 
reason to suspect the accuracy of the observations within this phase range. 
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Fic. 1.—Blue light curve of FU Arae. 
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Fic. 2.— Yellow light curve of FU Arae. 


The range in magnitude is very much greater than found by Miss Swope. 
At primary minimum it is 0™-97 in blue light and o™-8g in yellow light. At 
secondary minimum the respective values are o™-20 and o™-23. Thus primary 
minimum is much deeper than secondary. So far as can be seen from the light 
curve the latter occurs at 0-5 phase. 
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